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In this paper we investigate the origin of the mid-infrared (IR) hydrogen recombination lines for a sample of 
114 disks in different evolutionary stages (full, transitional and debris disks) collected from the Spitzer archive. 
We focus on the two brighter H I lines observed in the Spitzer spectra, the H 1(7-6) at 1 231 pm and the H 1(9-7) 
at 11.32/im. We detect the H 1(7-6) line in 46 objects, and the H 1(9-7) in 11. We compare these lines with the 
other most common gas line detected in Spitzer spectra, the [Ne II] at 12.81 pm. We argue that it is unlikely that 
the H I emission originates from the photoevaporating upper surface layers of the disk, as has been found for 
the [Ne II] lines toward low-accreting stars. Using the H l(9-7)/H 1(7-6) line ratios we find these gas lines are 
likely probing gas with hydrogen column densities of 10 10 -10 n cm -3 . The subsample of objects surrounded 
by full and transitional disks show a positive correlation between the accretion luminosity and the H I line 
luminosity. These two results suggest that the observed mid-IR H I lines trace gas accreting onto the star in 
the same way as other hydrogen recombination lines at shorter wavelengths. A pure chromospheric origin of 
these lines can be excluded for the vast majority of full and transitional disks. We report for the first time the 
detection of the H 1(7-6) line in eight young (< 20 Myr) debris disks. A pure chromospheric origin cannot be 
ruled out in these objects. If the H 1(7-6) line traces accretion in these older systems, as in the case of full and 
transitional disks, the strength of the emission implies accretion rates lower than l(r m MWyr. We discuss some 
advantages of extending accretion indicators to longer wavelengths, and the next steps required to pin down the 
origin of mid-infrared hydrogen lines. 

Subject headings: accretion, accretion disks — circumstellar matter — infrared: stars — line: identification 
— stars: activity 


1. INTRODUCTION 

The analysis and characterization of circumstellar disks, 
which provide the raw material to form planets, are crucial for 
our understanding of disk dispersal and planet formation pro¬ 
cesses. Observations and modeling of the dust content in pro¬ 
toplanetary disks are nowadays relatively abundant (see Testi 
et al. 2014 for a review). The gaseous component of circum¬ 
stellar disks is instead less well characterized (e.g., Dutrey 
et al. 2014 for a review) even if major advances have been 
made in the past few years thanks to space observatories such 
as Spitzer (e.g., Pontoppidan et al. 2014) and Herschel (e.g.. 
Dent et al. 2013). Gas plays a pivotal role in the structure 
and evolution of disks representing the bulk of the disk mass 
(e.g., Alexander et al. 2014). It is then necessary to extend 
our knowledge of the disk gas content in order to improve our 
understanding of the natal environments of planets. 

The sensitivity of the InfraRed Spectrograph (IRS) on the 
Spitzer Space Telescope led to the discovery of gas lines that 
can be used to characterize the immediate surroundings of 
young stellar objects (YSOs). For instance, the rich water 
spectrum of YSOs has been attributed to emission arising 
from warm (~600 K) layers of water vapor out to a few AUs 
from the central star (Carr & Najita 2008, 2011, Salyk et al. 
2008, Pontoppidan et al. 2010, Salyk et al. 2011, Banzatti et 
al. 2013). 

Among the wide forest of lines revealed by Spitzer the most 
common one is the singly ionized neon, [Ne II] at 12.81 pm 
(e.g., Pascucci et al. 2007, Lahuis et al 2007, Giidel et al. 

2010) . Different studies of this line show that objects with 
known jets/outflows have higher [Ne II] luminosity than ob¬ 
jects without (Giidel et al. 2010, Baldovin-Saavedra et al. 

2011) . Moreover, follow up observations of this line with 
ground-based high-resolution spectrographs have shown that 
[Ne II] is predominantly tracing an outflow in high-accreting 
stars, while its origin is from photoevaporative disk winds 
driven by high-energy photons from the central star in lower 
accretors (Pascucci & Sterzik 2009, van Boekel et al. 2009, 
Sacco et al. 2012, Baldovin-Saavedra et al. 2012). 

Other gas lines that are commonly used to characterize ob¬ 
jects undergoing accretion are hydrogen recombination lines. 

H I lines from the near-ultraviolet (UV) to the near-infrared 
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(IR) have been extensively used as indicators of the rate at 
which the matter from the inner circumstellar disk is accret¬ 
ing onto the star (e.g., Muzerolle et al. 2001, Calvet et al. 
2004, Natta et al. 2006, Alcala et al. 2014, Edwards et al. 
2013 among many others). Recently, Salyk et al. (2013) using 
NIRSPEC at the Keck II telescope and CRIRES at the Very 
Large Telescope have introduced the H I Pf/3 line at 4.65/im 
as accretion diagnostic. This is to date the longest wavelength 
hydrogen line proposed as accretion indicator. 

As early as 2007, Pascucci et al. (2007) reported also the 
detection of H I lines in the Spitzer spectrum of RX1852.3- 
3700. Subsequently, Ratzka et al. (2007) and Najita et al. 
(2010) reported hydrogen line detections in both the low- and 
high-resolution Spitzer spectrum of TWHya. The analysis of 
these two objects did not constrain the origin of these lines and 
left several possibilities open: stellar chromosphere, hot disk 
atmosphere, accretion shock regions, or internal wind shock 
very close to the stellar surface (Hollenbach & Gorti 2009). 
More recently, Carr & Najita (2011) reported the detection of 
the IR hydrogen lines in 10 T Tauri stars but did not discuss 
their origin. 

Here, we improve upon studies of gas lines from circum¬ 
stellar disks by analyzing hydrogen lines observed in the 
Spitzer/IRS spectra of a large sample of objects. In particular 
we focus on the strongest of the Humphreys Series lines in the 
mid-infrared, the Hu a at 1231 pm (hereafter H 1(7-6)) and the 
H 1(9-7) line at 11.32/rm. With the aim to investigate the ori¬ 
gin of these mid-IR hydrogen lines, we compare their proper¬ 
ties both with hydrogen lines observed at shorter wavelengths, 
mainly tracing accretion, and with the other most common gas 
line observed in the Spitzer spectra ([Ne II]) which is instead 
an outflow or disk wind tracer depending on the mass accre¬ 
tion rate. 

The paper is organized as follows. In Section 2 the source 
selection and data reduction are explained. In Section 3 we 
describe the method used to account for the water contamina¬ 
tion of the H I lines, and how we measure the H I and [Ne II] 
line fluxes and luminosities. In Section 4 we discuss the origin 
of the H I mid-IR lines, highlighting the interesting presence 
of H I lines in a sample of debris disks. The summary of the 
results and an outlook are given in Section 5. 

2. SAMPLE SELECTION AND DATA REDUCTION 
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Using the Spitzer Heritage Archive 1 we downloaded and 
re-reduced all IRS Short-High (SH) spectra of disks-bearing 
young stars acquired in the Standard Staring Mode and which 
have multiple on-source exposures (N-cycle >1) and sky ob¬ 
servations. The latter two restrictions were applied to better 
identify and correct for bad pixels as summarized below. A 
list of the 114 objects identified in this way and with a signal- 
to-noise on the continuum greater than 5 is reported in Ta- 

1 http ://sha. ipac. caltech. edu/applications/Spitzer/SHA/ 


ble[Q This table also contains the source coordinates, the As¬ 
tronomical Observation Request number (AOR), the Identifi¬ 
cation Number (ID) of the Spitzer program the object belongs 
to, and the exposure time per frame times the number of ex¬ 
posures. We also report in Table lATI the sources’ properties, 
and in Table I A2l the mid-IR line fluxes, errors and continuum 
levels around the three line of interest. 


TABLE 1 

All objects reduced: Observing log 


Name 

RA 

DEC 

AOR 

Program ID 

Time x Ncycles 

HD 377 

0h08m25.76s 

+6d37m00.3s 

13462272 

148 

120 x 3 

HD 12039 

Ih57m48.97s 

-21d54m05.2s 

13461760 

148 

120 x 4 

HD 17925 

2h52m32.13s 

-12d46m 11.0s 

9780480 

148 

30 x 5 

HD 19668 

3h09m42.28s 

-9d34m46.4s 

13462784 

148 

120 x 3 

LkHa 326 

3h30m44.01s 

+30d32m47.0s 

27063552 

50641 

30 x 12 

IC348-67 

3h43m44.62s 

+32d08ml7.9s 

22850816 

40247 

30 x 100 

IC348-31 

3h44ml8.16s 

+32d04m57.0s 

22849280 

40247 

30 x 16 

IC348-72 

3h44m22.58s 

+32d01m53.7s 

22851840 

40247 

30 x 100 

IC348-68 

3h44m28.51s 

+31d59m54.1s 

22851328 

40247 

30 x 100 

IC348-55 

3h44m31.36s 

+32d00ml4.7s 

22850304 

40247 

30 x 80 

IC348-2 

3h44m35.36s 

+32dl0m04.6s 

22847744 

40247 

30 x 4 

IC348-6 

3h44m36.94s 

+32d06ra45.4s 

22848256 

40247 

30 x 16 

IC348-37 

3h44m37.99s 

+32d03m29.8s 

22849792 

40247 

30 x 30 

IC348-133 

3h44m41.74s 

+32dl2m02.4s 

22852352 

40247 

30 x 100 

IC348-21 

3h44m56.15s 

+32d09ml5.5s 

22848768 

40247 

30 x 30 

LkHa 330 

3h45m48.28s 

+32d24ml 1.9s 

27063040 

50641 

30 x 7 

HD 25457 

4h02m36.74s 

-0dl6m08.1s 

9779712 

148 

6x4 

FM Tau 

4hl4ml3.58s 

+28dl2m49.3s 

24403968 

30300 

120 x 2 

FN Tau 

4hl4ml4.59s 

+28d27m58.0s 

27062528 

50641 

30 x 6 

CW Tau 

4hl4ml7.00s 

+28dl0m57.8s 

27062272 

50641 

30 x 6 

FP Tau 

4hl4m47.30s 

+26d46m26.4s 

26321408 

50498 

30 x 8 

CX Tau 

4hl4m47.86s 

+26d48m 11.0s 

23551744 

40338 

30 x 21 

CYTau 

4hl7m33.72s 

+28d20m46.8s 

27061504 

50641 

30 x 14 

DD Tau 

4hl8m31.13s 

+28dl6m28.9s 

26320384 

50498 

30 x 2 

BP Tau 

4hl9ml5.84s 

+29d06m27.0s 

14548224 

20363 

30 x 12 

DE Tau 

4h21m55.67s 

+27d55m06.2s 

27060992 

50641 

30 x 10 

FT Tau 

4h23m39.18s 

+24d56ml4.3s 

27060224 

50641 

30 x 14 

04216+2603 

4h24m44.58s 

+26dl0ml4.2s 

27059968 

50641 

30 x 12 

IP Tau 

4h24m57.08s 

+27dllm56.3s 

24404224 

30300 

120 x 5 

DG Tau 

4h27m04.70s 

+26d06m 15.8s 

14547968 

20363 

6 x 12 

DHTau 

4h29m41.57s 

+26d32m58.2s 

26320128 

50498 

30 x 10 

UX TauA 

4h30m04.00s 

+18dl3m49.4s 

18013952 

30300 

30 x 30 

DKTau 

4h30m44.28s 

+26d01m24.6s 

14548736 

20363 

30 x 8 

HK Tau 

4h31m50.58s 

+24d24ml7.8s 

27059200 

50641 

30 x 14 

FZ Tau 

4h32m31.77 s 

+24d20m02.6s 

27058688 

50641 

30 x 6 

HN Tau 

4h33m39.36s 

+17d51m52.3s 

27058176 

50641 

30 x 6 

DL Tau 

4h33m39.08s 

+25d20m38.1s 

27058432 

50641 

30 x 6 

DM Tau 

4h33m48.72s 

+18d 10m 10.0s 

18014720 

30300 

30 x 30 

AA Tau 

4h34m55.44s 

+24d28m53.4s 

14551552 

20363 

30 x 12 

DN Tau 

4h35m27.37s 

+24dl4m58.9s 

24404480 

30300 

30 x 9 

DO Tau 

4h38m28.58s 

+26dl0m49.4s 

14548480 

20363 

30 x 8 

LkCa 15 

4h39ml7.80s 

+22d21m03.5s 

23551488 

40338 

30 x 21 

04385+2550 

4h41m38.82s 

+25d56m27.0s 

27057408 

50641 

30 x 6 

DP Tau 

4h42m37.69s 

+25dl5m37.3s 

27057152 

50641 

30 x 6 

DQ Tau 

4h46m53.06s 

+17d00m00.4s 

26321152 

50498 

30 x 2 

DR Tau 

4h47m06.22s 

+16d58m42.9s 

27067136 

50641 

6 x 12 

DS Tau 

4h47m48.60s 

+29d25m 11.6s 

26320896 

50498 

30 x 4 

UY Aur 

4h51m47.38s 

+30d47ml3.9s 

14551040 

20363 

6 x 10 

St 34 

4h54m23.68s 

+17d09m53.5s 

18016256 

30300 

30 x 40 

GM Aur 

4h55ml0.98s 

+30d21m59.5s 

18015488 

30300 

30 x 10 

SUAur 

4h55m59.38s 

+30d34m01.5s 

27066880 

50641 

6 x 10 

HD 32297 

5h02m27.44s 

+7d27m39.7s 

25677568 

50150 

30 x 5 

V836 Tau 

5h03m06.60s 

+25d23ml9.7s 

23552000 

40338 

30 x 21 

RW Aur 

5h07m49.56s 

+30d24m05.0s 

14547200 

20363 

6 x 13 

HD 35850 

5h27m04.76s 

-1 Id54m03.5s 

9777920 

148 

30 x 6 

HD 37484 

5h37m39.62s 

-28d37m34.7s 

9780224 

148 

30 x 10 

NGC2068/ 2071-H9 

5h46ml 1.86s 

+0d32m25.9s 

22853888 

40247 

30 x 100 

NGC2068/ 2071-SI 

5h46m 18.89s 

-0d05m38.2s 

22854912 

40247 

30 x 50 

NGC2068/ 2071-HI3 

5h46m22.44s 

-0d08m52.6s 

22854400 

40247 

30 x 80 

AO Men 

6hl8m28.24s 

-72d02m41.7s 

5458688 

148 

120 x 2 

HD 61005 

7h35m47.47s 

-32dl2ml4.2s 

13462528 

148 

120 x 3 
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TABLE 1 — Continued 


Name 

RA 

DEC 

AOR 

Program ID Time x Ncycles 

HD 69830 

8hl8m23.95s 

-12d37m55.8s 

12710656 

41 

30 x 2 

SX Cha 

10h55m59.76s 

-77d24m40. Is 

27066624 

50641 

30 x 7 

TW Cha 

10h56m30.52s 

-77dllm39.4s 

27066368 

50641 

30 x 28 

SZCha 

10h58ml6.77s 

-77dl7ml7.1s 

22846208 

40247 

30 x 4 

TW Hya 

1 Ih01m51.91s 

-34d42ml7.0s 

18017792 

30300 

30 x 2 

CS Cha 

1 Ih02m24.91s 

-77d33m35.7s 

18021632 

30300 

30 x 50 

Cedi 10-IRS2 

1 Ih06ml5.41s 

-77d21m56.8s 

18022400 

30300 

30 x 10 

Cedi 10-IRS4 

1 Ih06m46.58s 

-77d22m32.6s 

18019328 

30300 

30 x 36 

Sz 18 

Ilh07ml9.15s 

-76d03m04.8s 

18020864 

30300 

30 x 100 

VW Cha 

Ilh08m01.47s 

-77d42m28.7s 

27066112 

50641 

30 x 6 

Sz 27 

1 Ih08m39.05s 

-77dl6m04.2s 

22847232 

40247 

30 x 100 

VZ Cha 

1 Ih09m23.79s 

-76d23m20.8s 

27065856 

50641 

30 x 8 

WX Cha 

1 Ih09m58.74s 

-77d37m09.0s 

27065600 

50641 

30 x 12 

Hen 3-600A 

1 Ihl0m27.53s 

-37d31m53.6s 

18018560 

30300 

30 x 2 

XX Cha 

1 lhl lm39.66s 

-76d20ml5.3s 

27065344 

50641 

30 x 14 

RX Jill 1.7-7620 

llhllm46.31s 

-76d20m09. Is 

5451776 

148 

30 x 6 

CHX22 

1 Ihl2m42.69s 

-77d22m23.1s 

22846720 

40247 

30 x 60 

1RXS J121236.4-552037 

12hl2m35.77s 

-55d20m27.4s 

5460992 

148 

120 x 6 

1RXS J122233.4-533347 

12h22m33.23s 

-53d33m48.9s 

13463296 

148 

120 x 6 

TWA 10 

12h35m04.30s 

-41d36m39.0s 

22934528 

40338 

30 x 10 

Sz 50 

13h00m55.38s 

-77dl0m22.1s 

27065088 

50641 

30 x 14 

1RXS J130153.7-530446 

13h01m50.69s 

-53d04m58.1s 

13462016 

148 

120 x 6 

1RXS J132207.2-693812 (MPMus) 

13h22m07.54s 

-69d38ml2.1s 

5451264 

148 

6x3 

1RXS J133758.0-413448 

13h37m57.30s 

-41d34m42.0s 

5460224 

148 

120 x 6 

HD 119269 

13h43m28.54s 

-54d36m43.5s 

5461248 

148 

120 x 2 

HD 134319 

15h05m49.91s 

+64d02m50.0s 

9779968 

148 

30 x 6 

V343Nor 

15h38m57.56s 

-57d42m27.3s 

5458944 

148 

120 x 2 

HTLup 

15h45ml2.87s 

-34dl7m30.6s 

27064832 

50641 

6 x 10 

GQ Lup 

15h49ml2.11s 

-35d39m04.8s 

27064576 

50641 

30 x 8 

HD 141943 

15h53m27.30s 

-42dl6m00.8s 

5462528 

148 

120 x 4 

HD 142361 

15h54m59.86s 

-23d47ml8.2s 

5461504 

148 

120 x 12 

IM Lup 

15h56m09.21s 

-37d56m05.9s 

27064320 

50641 

30 x 6 

RULup 

15h56m42.31s 

-37d49ml5.5s 

27064064 

50641 

30 x 5 

HD 143006 

15h58m36.93s 

-22d57ml5.4s 

9777152 

148 

6x8 

RX J1600.6-2159 

16h00m40.57s 

-22d00m32.3s 

17145344 

148 

120 x 9 

EX Lup 

16h03m05.49s 

-40dl8m25.3s 

27063808 

50641 

30 x 6 

RX J1612.6-1859 

16hl2m40.51s 

-18d59m28.3s 

5457664 

148 

120 x 4 

RXJS J161410.6-230542 

16hl4ml 1.08s 

-23d05m36.2s 

5453824 

148 

30 x 2 

1RXS J161458.4-275013 

16hl4m59.17s 

-27d50m23.0s 

17144832 

148 

120 x 9 

DoAr 21 

16h26m03.03s 

-24d23m36.4s 

24403456 

30300 

30 x 6 

DoAr 24E 

16h26m23.42s 

-24d21m00.5s 

27062784 

50641 

30 x 5 

DoAr 25 

16h26m23.69s 

-24d43m 14.0s 

24403712 

30300 

120 x 4 

SR 21 

16h27m 10.20s 

-24d 19m 16.0s 

24403200 

30300 

6 x 10 

Haro 1-16 

16h31m33.46s 

-24d27m37.4s 

27062016 

50641 

30 x 7 

RNO 90 

16h34m09.17s 

-15d48m 16.8s 

27061760 

50641 

6x9 

WaOph6 

16h48m45.63s 

-14dl6m36.0s 

27060736 

50641 

30 x 6 

EC 82 

18h29m56.90s 

+Idl4m46.5s 

27059712 

50641 

30 x 8 

RX J1842.9-3532 

18h42m57.97s 

-35d32m42.7s 

5451521 

148 

30 x 5 

RXJ 1852.3-3700 

18h52ml7.30s 

-37dOOm 11.9s 

5452033 

148 

120 x 2 

AS 353A 

19h20m30.98s 

+ 1 Id01m54.8s 

14547456 

20363 

30 x 8 

V1331 Cyg 

21h01m09.20s 

+50d21m44.5s 

14547712 

20363 

30 x 8 

HD 202917 

21h20m49.96s 

-53d02m03.2s 

9778176 

148 

30 x 5 

HD 216803 

22h56m24.05s 

-31d33m56.1s 

9777664 

148 

30 x 5 


The IRS SH module uses a cross-dispersed echelle design 
and covers the spectral range 9.9-19.6 /im with a nominal re¬ 
solving power of R^600 (Houck et al. 2004). The data re¬ 
duction was performed starting from the droopres products 
following the procedure outlined in Pascucci et al. (2013). 
In brief, we computed the average of the on-source and sky 
exposures per nod position and rejected bad pixels using the 
noise statistics from multiple cycles. We then created back¬ 
ground subtracted images and corrected for the additional 
known bad pixels reported in the IRS BCD mask. The one 
dimensional spectrum was extracted from the background- 
subtracted pixel-corrected images using the full aperture ex¬ 
traction routine in SMART (Higdon et al. 2004). We applied 
the same procedure to our sources and to the Spitzer flux cal¬ 
ibrator £ Dra. From the reduction of 10 different exposures 
of £ Dra and the known stellar model atmosphere we created 


two one-dimensional spectral response functions (one at each 
nod position). Finally, we divided the source spectra by the 
spectral response function. A compilation of four Spitzer/ IRS 
SH spectra analyzed in this paper is shown in Figure [Q 
The calibrated spectra of 54 sources belonging to 
the program IDs 40247, 20363, 148 and 41 can be 
downloaded from the IRSA Spitzer Science Center at 
http://irsa.ipac.caltech.edu/data/SPITZER/Disks_SH_spectra 
The total sample includes objects classified according to 
their general disk properties. We use the generic term Full 
Disks (F) for those disks with strong excess emission above 
photospheric levels at infrared wavelengths, starting in the 
near-IR. Transitional Disks (T) are broadly defined as disks 
with reduced near- and mid-IR excess fluxes but excess emis¬ 
sion comparable to full disks beyond ~20//m, indicative of 
a depleted inner region. Finally, Debris Disks (D) are de- 
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fined in this paper as having pure photosphere emission out 
to about 10 pm but can have relatively small excess emission 
in the long portion of the mid-IR spectrum. This disk classifi¬ 
cation collected from the literature is based on the dust prop¬ 
erties and the continuum spectral energy distributions, and it 
is mainly taken from Furlan et al. (2006, 2011), Espaillat et 
al. (2012), Silverstone et al. (2006), Chen et al. (2014), as 
reported in Table IA1I Most of the debris disks included in 
our study belong to the recent catalog published by Chen et 
al. (2014). Full and transitional disks are primordial, gas-rich 
disks, while debris disks are second-generation disks, domi¬ 
nated by dust. In total our sample includes 68 full, 21 transi¬ 
tional, and 25 debris disks. In Figure|2] we show the distribu¬ 
tion in mass, luminosity, effective temperature and accretion 
luminosity (see Sect. l3.2l > for the total sample, and for the ob¬ 
jects where the H 1(7-6) line is detected. A K-S two sample 
test of the stellar properties shows the two dataset have most 
likely the same parent distribution. 



Fig. 1.— Spitzer /IRS SH spectra of four disks. For plotting purposes 
offsets have been applied as follows: -0.1 Jy at VWCha (water-contaminated 
full disk), -0.6 Jy at AS 353A (water-free full disk), +0.3 Jy at GM Aur 
(transitional disk) and 0.0 Jy at HD 35850 (debris disk). The vertical dashed 
lines show the location of the three lines mainly used through this paper: 
H 1(7-6) at 12.37^m, H 1(9-7) at 11.32/im and [Ne II] at 12.81^m. 


3. ANALYSIS AND IMMEDIATE RESULTS 

In the following, we summarize the main properties of our 
sample (Sect. 13.1k compute accretion luminosities for each 
object in a self consistent manner (Sect. l3T2k and describe in 
detail how we calculate [Ne II] and H I fluxes taking into ac¬ 
count any possible contamination from water lines (Sects. [3731 
and 13.4k The last two sections (Sects. 1331 and [331 discuss 
detections and correlations between line fluxes and adjacent 
continuum emission. 

3.1. Subsample ofH I and [Ne II7 sources: star and disk 
properties 

The properties of the 114 objects are summarized in Ta¬ 
ble IA1I For each object we provide the stellar association 
it belongs to, the distance, spectral type (SpTy), visual ex¬ 
tinction (Ay) as well as the effective temperature (T e ff) ob- 




l°g L,l.r (Mg) 




Fig. 2.— Distribution of mass, luminosity, effective temperature and accre¬ 
tion luminosity for all the objects in our sample (gray), and only the objects 
where H 1(7-6) line has been detected (blue). 


tained using the T e ff-SpTy correlation of Kenyon & Hartmann 
(1995) and Fuhman et al. (2003). Ha equivalent widths and 
stellar luminosities (L slar ) have been collected from the litera¬ 
ture using the references reported in Table lAll while the mass 
(M, fflr ), and radius (R srar ) were derived from the HR diagram 
and the Siess et al. (2000) evolutionary tracks. 

For multiple star systems the name we report in Table lAll is 
that of the primary star, which is the one with the earliest spec¬ 
tral type. However, in the table caption we list the separation 
of any other stellar companion within 10" of the primary star. 
Given that the IRS/SH slit width is 4.7", there are 15 objects 
with companions closer than 2.5" that can contaminate the ex¬ 
tracted IRS spectra. However, McCabe et al. (2006) show that 
only in 25% of the binary systems the secondary flux becomes 
larger than the primary flux at wavelengths >2 pm. Based on 
this, we assign the Spitzer spectrum always to the primary star 
in the 15 unresolved multiple systems. 

We emphasize that the sample collected is biased in several 
ways: i) the number of objects with sky observations and mul¬ 
tiple on-source exposures (see Sect. [2]), ii) we only consider 
objects where the continuum has been detected and the signal- 
to-noise ratio of continuum level is higher than 5. Hence, we 
do not attempt any statistical study based on this sample. 


3.2. Accretion Properties 

To test whether mid-IR H I lines trace the accretion shock 
region and/or the accretion flow, we compare H I line lumi¬ 
nosities with stellar accretion luminosities (L acc ). One of the 
most direct methods to measure L acc is through modeling of 
the UV excess continuum emission produced by the accretion 
shock (e.g. Gullbring et al. 1998). However, only 17 objects 
in our sample have such measurements (AA Tau, CS Cha, 
DE Tau, DK Tau, DN Tau, DM Tau, DR Tau, GM Aur, 
HN Tau, FkCa 15, RX 132207.2-693812, TW Hya, SZ Cha, 
Sz 27, FkHa 330, RX J1852.3-3700 and RX J1842.9-3532, 
see Ingleby et al. 2013 and Manara et al. 2014). To measure 
accretion luminosities in a consistent way for all the objects 
in our sample we will use the Ha line luminosity. Even if this 
method allows an homogeneous estimate of the accretion lu- 
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minosity, we stress that the hydrogen lines in general, and the 
Ha line in particular, have several contributing mechanisms: 
accretion, stellar winds, jets, chromosphere (e.g. Hartmann et 
al. 1990, Calvet & Hartmann 1992). While this line is known 
to be an imperfect tracer of accretion in young stars, Rigliaco 
et al. (2012), among others, have shown that L acc obtained 
from the UV excess emission and via indirect methods (e.g., 
Ha line luminosity) agree to within an order of magnitude. 

To measure the accretion luminosity we first compute the 
Ha luminosity from the Ha equivalent width, the R- and/or 
/-band magnitude, and the extinction (Ay) for each object 
(e.g., Hartigan et al. 1995). The data were dereddened us¬ 
ing Mathis’s reddening law for R v =3.1 (Mathis 1990). Next, 
we convert the Ha luminosity into L acc using the relationship 
found by Rigliaco et al. (2012). Although this method pro¬ 
vides self-consistent values for L acc we are aware that the not- 
simultaneous observation of Ha, photometry, and Ay might 
introduce additional uncertainties in our estimates of L acc . In 
Sect. 14.21 we describe how we account for these sources of 
uncertainty. 

3.3. [Ne II/ Line Fluxes and Luminosities 

Water is largely detected in protoplanetary disks (e.g., Pon- 
toppidan et al. 2010, Riviere-Marichalar et al. 2012), and 
its emission can contaminate other emission lines, as shown 
in the following sections. This is not the case for the [Ne II] 
line, that is quite isolated from other HjO transitions. For this 
water-isolated and spectrally unresolved gas line we compute 
the flux by fitting a Gaussian profile to the observed line. We 
define the continuum level around the line as the mean value 
obtained sampling a 0.1pm wide region shortward and long- 
ward of the transition, and the uncertainty on the measured 
flux is the rms in these regions. The line flux and the associ¬ 
ated la uncertainty are reported in Table lA2l for each detec¬ 
tion. In the same table we provide 3er upper limits to the flux 
when [Ne II] lines are not detected. Line luminosities are then 
computed using the distances reported in Table lAll 

We have co mpare d the measured [Ne II] line fluxes (as re¬ 
ported in Table IA2I) with those reported by other authors us¬ 
ing either different data reduction techniques (Pascucci et al. 
2007, Lahuis et al. 2007, Espaillat et al. 2013), and/or includ¬ 
ing any possible unresolved contribution from other species 
(HjO and OH, Carr & Najita 2011). We find that the mea¬ 
sured line fluxes and upper limits are in agreement within the 
errors (see Fig.0. 

3.4. H I Line Fluxes and Luminosities 

The method used to compute the [Ne II] line fluxes and lu¬ 
minosities could be in principle applied to the H I lines. How¬ 
ever, because the H 1(7-6) and H 1(9-7) transitions are very 
close to water lines we must subtract any possible contami¬ 
nation from water before measuring line fluxes. We proceed 
as follows. In Sect. 13.4.11 we identify the objects that lack 
detectable water emission (hereafter, ’water-free’ disks 2 ) and 
the objects where water is detected in their spectra (’water- 
contaminated’ disks). For the water-free disks we compute 
line fluxes as in Sect. 13.31 For the water-contaminated disks 
we discuss a procedure to subtract water emission from the 
H I lines (Sect. 13.4.2b and then we measure the residual H I 
fluxes after water subtraction. 

2 We note that these objects might possibly show water contribution in 
higher S/N or higher resolution spectra. 



Fig. 3.— Comparison between the [Ne II] line fluxes measured fol¬ 
lowing the data reduction chain explained in Sect. [2] and the line 
fluxes obtained from the literature using different data reduction 
techniques. 

3.4.1. Water-contaminated vs water-free disks 

To search for water emission in the source spectra we used 
the two water line complexes at 15.17pm and 17.22pm iden¬ 
tified by Pontoppidan et al. (2010). These lines were selected 
because they are strong and more isolated than most other 
HaO lines. Following Pontoppidan et al. (2010), we define 
as water contaminated those spectra in which both the H 2 O 
lines are detected above 3a. We do not find any case where 
one of the two lines is detected above 3a, and the other one is 
not detected. 

Among the 114 disks in our sample, 33 are water- 
contaminated. The remaining objects have water-free spec¬ 
tra. We show in Figure [4] a water-contaminated spectrum 
(DR Tau), and a water-free spectrum (DE Tau) in the region 
around the 15.17pm and 17.22pm water lines. For compari¬ 
son we also show a synthetic water spectrum modeled assum¬ 
ing a gas temperature of 600 K and a water column density 
of 10 18 cm -2 (see next Section). Line fluxes, together with 
the uncertainties and continuum level next to the lines are re¬ 
ported in Table lA2l 

3.4.2. Correcting for water emission 

Water emission lines have been detected towards many pro¬ 
toplanetary disks using the Spitzer/IRS spectrograph (Pontop¬ 
pidan et al. 2014 for a review). Previous studies have shown 
that a simple model, which assumes a plane-parallel slab ge¬ 
ometry and gas in local thermal equilibrium (LTE), can re¬ 
produce relatively well the fluxes of the observed H 2 O lines 
(e.g., Carr & Najita 2011). However, we are aware that these 
assumptions ignore the non-LTE effects and the complex ge¬ 
ometric structure of disks (see Meijerink et al. 2009). We 
use this simple model as developed in Pascucci et al. (2013) 
to subtract any water contamination from the H I lines. The 
specific LTE model takes the molecular parameters from the 
HITRAN 2008 database (Rothman et al. 2009) and requires 
only 4 inputs: the gas temperature, the water column density, 
the radius of the emitting area, and the line width (which we 
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Fig. 4.— Examples of water contaminated (DR Tau), water-free 
(DE Tau) and synthetic H 2 O spectrum in the wavelength region 
around the two isolated H 2 O lines at 15.17 and 17.22 /tm. The tar¬ 
get spectra are shown in black, and the synthetic water spectrum is 
shown in red. 

assume to be the thermal line width, i.e., variable with tem¬ 
perature). We created a grid of synthetic spectra using five 
different gas temperatures (T=[400, 500, 600, 700, 800]) and 
water column densities (N ra /=[10 i7 , 5xl0 17 , 10 18 , 5xl0 18 , 
10 19 cm -2 ]) and assumed the same projected emitting area of 
radius 1 AU. The chosen inputs span the range of parameters 
inferred for the warm disk atmosphere of T Tauri stars (Carr 
& Najita 2011, Salyk et al. 2011). 

The water lines contaminating the hydrogen lines are cen¬ 
tered at ~12.39/rm and ~11.31/xm for the H 1(7-6) and 
H 1(9-7) transition, respectively. In the following we will re¬ 
fer to these lines as ‘T^O-contaminating”. To subtract off 
the water contribution to the observed H I lines we have 
searched the synthetic spectra for isolated EDO lines with 
similar properties as the H 2 0-contaminating ones. The two 
water lines that we have identified will be called in the fol¬ 
lowing “EDO-isolated” lines and have the following charac¬ 
teristics: they have similar lower energies (E i ow ) and Ein¬ 
stein A-coefficients (A,,/) (within 15% at most) as the EEO- 
contaminating lines; EDO-isolated line is within 1 //m from 
the corresponding contaminating line, so that the continuum 
level does not substantially change. Following these crite¬ 
ria we find that the EDO-isolated line at 13.13/rm is the best 
match for the EDO-contaminating line at 12.39/rm, while the 
EDO-isolated at 11.96/im has very similar properties as the 
EDO-contaminating at 11.31 //m (see Table |2}. We show in 
Fig. [5]the ratio between the H 2 O contaminating/isolated line 
fluxes for all the temperature and column density of the syn¬ 
thetic spectra. The ratios are quite constant among all the syn¬ 
thetic spectra, with EDO(12.39/im)/EDO(13.13/im)~3.5 and 
EDO(lD32/rm)/EDO(l 1.97/xm)~1.0, and vary within a factor 
1.5 at most over our grid of synthetic spectra. 

Once we identified the best water proxies, we searched all 
the spectra for EDO-isolated lines. In the SH-Spitzer/IRS 
spectra none of the 83 water-free spectra shows these lines, 
as expected because the H 2 O lines become weaker at shorter 
wavelengths (e.g., Salyk et al. 2011). In the 33 water- 
contaminated disks, we subtracted the water contribution 
from the H I line fluxes as follows: if the EDO-isolated 
line is detected above the 3er level, we use it and its ra¬ 
tio to the nearby EDO-isolated line to compute the EDO- 


TABLE 2 

FLO Line Properties 


Line 

Wavelength 

Aul 

E, ow 


(/im) 

(s -1 ) 

(cm- 1 ) 

H 2 O-contaminating 

11.3241 

1.576 

2746.02 

H 2 O-contaminating 

12.3962 

7.665 

3211.21 

H20-isolated 

11.9681 

1.281 

2629.33 

H20-isolated 

13.1325 

6.788 

3244.60 


contaminating flux (as EDO(12.39/rm)=3.5xEDO(13.13/Tm) 
and EDO(l 1.32/im)=EDO(l 1.97/rm), respectively for H 1(7- 
6) and H 1(9-7)), and then subtract it from the measured HI 
line flux. If the H 2 0-isolated is not detected, we subtract its 
3a upper limit from the observed H I flux. After subtracting 
the water contamination from the observed line flux, we ob¬ 
tain the actual hydrogen line flux. The final line fluxes after 
water subtraction and errors are reported in Table lA2l We also 
show in parenthesis the total flux observed before subtracting 
the water contribution in the objects where the H I lines are 
detected. In the case of H I non-detections we provide the 
3a upper limit to the flux. Any other water contribution, if 
present, is within the 10% absolute flux uncertainty. As for 
the [Ne II], we have compared our H 1(7-6) line flux measure¬ 
ments with the few already reported in the literature (Pascucci 
et al. 2007, Carr & Najita 2011). We find that the reported val¬ 
ues are in agreement within the errors and the absolute flux 
uncertainty. 


o 

e” 

O 
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Fig. 5.— H 2 O contaminating/isolated line ratios: 
crosses for H 2 0(11.32/tm)/H20(l 1.97/rm) and squares for 
H 2 0(12.39/2m)/H20(13.13/rm). The ratios are computed for all the 
temperatures and column densities explored with the model grid. At 
each temperature there are 5 ratios, one for each value of the column 
density. 


3.5. Line Detections Summary 

We detect the [Ne II] line at 12.81 /im in 49 objects in our 
sample. Among the 67 objects where the [Ne II ] line is not 
detected, 37 are full disks, 4 are transitional disks and 26 are 
debris disks. One debris disk has a [Ne II] line detection 
(Rigliaco et al. in preparation). After correcting for water 
contamination, we identify 46 objects with an H 1(7-6) detec¬ 
tion above 3a, 11 of them also show the weaker H 1(9-7) line. 
In 6 spectra we detect all three lines of interest: the [Ne II], the 
H 1(7-6), and the H 1(9-7) line. Line detections are summa- 
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TABLE 3 
Line Detections 


Detected Lines 

Total number 

Relative % 


of objects 

F 

T 

D 

[Ne II] 

49 

67 

31 

2 

H 1(7-6) 

46 

57 

26 

17 

H1(9-7) 

11 

55 

27 

18 

[Ne II]+H 1(7-6) 

24 

50 

46 

4 

H I(7-6)+H 1(9-7) 

11 

55 

27 

18 

[Ne n]+H I(7-6)+H 1(9-7) 

6 

50 

50 

0 


Note. — The second column shows the total number of 
objects where the lines listed in Column 1 have been detected. 
The 3rd, 4th and 5th columns report the relative percent of 
these objects in full, transitional and debris disks respectively. 


rized in Table[3] We note that all objects where water emission 
is detected are full disks and none of the transitional or debris 
disk spectra have detectable water emission. Thus, expand¬ 
ing upon previous results (e.g., Pontoppidan et al. 2010), our 
analysis shows that the absence of water lines in the SH mod¬ 
ule is typical to more evolved transitional and debris disks, 
while full disks can either have a water rich spectrum or not. 

Interestingly, we detect H 1(7-6) in 8 young (<20 Myrs) de¬ 
bris disks. We will discuss these objects in Sect. 14.51 

3.6. H I lines versus continuum luminosity 

Once we accounted for the water contribution and mea¬ 
sured the H 1(7-6) and H 1(9-7) line luminosities, we inves¬ 
tigated their behavior as a function of the continuum luminos¬ 
ity around the lines (Fig. [6]). As also seen in other IR lines 
(H 2 S(2), [Fe I] at 24/im, [Ne II], e.g., Lahuis et al. 2007), 
there is a correlation between the H I mid-IR lines and the 
continuum flux level next to the lines. This correlation could 
be influenced, but it is not exclusively driven, by the more 
favorable contrast between line and continuum luminosities 
where the continuum is lower (as demonstrated by finding the 
same increasing trend between the line to continuum vs con¬ 
tinuum flux). We have also investigated this correlation for 
different disk morphologies and find that debris disks have 
lower H I line luminosities and continuum luminosities than 
transitional and full disks. In order to understand whether 
these differences are statistically significant we have mea¬ 
sured the mean and standard deviation of line and continuum 
luminosities for the three different disk morphologies (shown 
as solid lines in Fig. [6]». From this analysis we conclude that, 
even if transitional disks appear having lower H I and contin¬ 
uum luminosities than full disks, their difference is within lcr 
uncertainty, and thus not statistically significant. On the other 
hand, debris disks show consistently lower luminosities, both 
in the continuum and in the lines. 

4. DISCUSSION 

4.1. Mid-IR H I lines: disk tracers? 

In the following we compare the [Ne II] line fluxes obtained 
from our reduction of the Spitzer spectra with those of the 
H I lines to test if these two lines could trace the same re¬ 
gion in/around the disk. In particular, we will test if [Ne II] 
and H I lines may trace the disk surface. As mentioned in 
Sect.[Q [Ne II] lines have been studied both with Spitzer and 
with ground-based spectrographs at higher spectral resolu¬ 
tion. These studies have shown that the [Ne II] line mostly 
traces jets/outflows in high-accreting stars and the hot disk 



Fig. 6. — H 1(7-6), H 1(9-7) and [Ne II] line luminosity versus 
continuum luminosity around the lines. Black-diamonds represent 
objects surrounded by a full disk, red-triangles are for transitional 
disks, blue-squares for debris disks. Downward pointing arrows are 
upper limits. In the left panel solid lines represent the mean value 
of line and continuum luminosities with the la uncertainty for the 
objects where both line and continuum are detected. In all cases, the 
mean value for the x-axis has been shifted by +1.5 for clarity. In 
the middle panel we show the mean in the line luminosity as dashed 
lines, given the large fraction of non-detections. We note again that 
the measure of the means is only based on the detections, and for 
the H 1(9-7) lines the upper limits are more numerous than the detec¬ 
tions. 


surface in lower accretors ( M acc < 5 x 10 -8 MQ/yr, e.g. Gtidel 
et al. 2010, Sacco et al. 2012). Unfortunately, similar ground- 
based followup studies have not yet been carried out for the 
mid-IR H I lines, hence a direct comparison of high-resolution 
spectra cannot be done. The vast majority of the objects in our 
sample (namely, 85%) for which we have computed L acc (and 
M acc ) can be classified as “low-accretors” according to the 
threshold value reported above (Giidel et al. 2010). Hence, 
for most of our sources a detection of [Ne II] in the Spitzer 
spectrum is most likely due to emission from the hot disk sur¬ 
face. 

The left panel of Figure^ shows a log-log plot of the H 1(7- 
6) line luminosity at 1 231pm versus the [Ne II] line luminos¬ 
ity at 12.81 pm . The luminosities of the two lines, when they 
are both detected, appear positively correlated, with Pearson 
correlation coefficient of ~0.9, and p-value~ 1 Q -6 . The cor¬ 
relation becomes weaker when upper limits are considered. 
However, this correlation could be affected by the underlying 
correlation between the line luminosities and the continuum 
luminosities around the line (see Figure [6] and discussion in 
Sect l3.6l) . To account for this possible effect we plot in Fig. [7] 
(right panel) the equivalent widths of the two observed lines. 
The [Ne II] and H 1(7-6) equivalent widths do not show any 
obvious correlation, and the Pearson correlation coefficient 
drops to ~0.4 (weak correlation) and p-value^0.04. This re¬ 
sult suggests that the positive relation between the H 1(7-6) 
and [Ne II] line luminosity is caused by the correlation be¬ 
tween the line and nearby continuum. 

To further test this finding, we compare the observed H I 
line luminosities with those predicted by models of disks ir¬ 
radiated by high-energy photons (Ercolano et al. 2008, Hol- 
lenbach & Gorti 2009, Ercolano & Owen 2010). In an EUV 
(13.6eV< hv <0.1 keV) irradiated disk the predicted Lhi( 7 - 6 )/ 
U[Neii] ratios are as small as ~0.008 (Hollenbach & Gorti 
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Fig. 7.— Top: H 1(7-6) versus [Ne II] line luminosity. We also show 
the prediction of these line luminosities for different theoretical mod¬ 
els. The black-dashed line (labeled as HG) shows the H I(7-6)/[Ne II] 
line relation for a disk irradiated by EUV photons (Hollenbach & 
Gorti 2009). The gray stars refer to X-ray irradiated primordial disk 
with logLx=29.3 (EOi), logL;c=30.3 (EO 2 ), and transitional disks il¬ 
luminated by logLx=30.3 erg/s with inner hole 8.3 AU (EO 3 ), 14.2 
AU (EO 4 ) and 30.5 AU (EO 5 ), as measured by Ercolano & Owen 
(2010). Bottom: comparison between the equivalent widths of the 
two lines. 

2009) because the radiative recombination rate of hydrogen 
is smaller than the electronic collisional excitation rate of 
[Ne II]. The observed line ratios are instead almost two orders 
of magnitude higher than the predicted ones (see Table lA2l >. 
A disk irradiated also by X-rays ( hv >0.1keV) still produces 
an H 1(7-6) line luminosity more than one order of magni¬ 
tude smaller than the ones observed (Ercolano et al. 2008, Er¬ 
colano & Owen 2010). The comparison between the observed 
line luminosity and different theoretical models shows that it 
is unlikely that the observed HI(7-6) lines trace the same hot 
disk atmosphere that produces the [Ne II] lines. 

In summary, both the lack of correlation between the [Ne II ] 
and the H 1(7-6) equivalent widths, and the comparison with 
model predictions suggest that the H 1(7-6) line is not tracing 
the same disk region as the [Ne II] line, and that the H I lines 
are unlikely disk gas tracers. 

4.2. Mid-IR H I lines: accretion indicators? 

Hydrogen recombination lines from the Balmer to the 
Paschen and Brackett series (in the UV, optical and near-IR 
wavelengths) have been extensively studied over the years, 
and have been convincingly found to trace infalling material 
onto the star, as demonstrated by their broad line profiles con¬ 
sistent with gas at almost free-fall velocity onto the star (e.g., 
Alcala et al. 2014). Nevertheless, these line profiles are also 
known to have other contributing mechanisms: stellar wind, 
jet, chromosphere. Recently, Salyk et al. (2013), using high- 
resolution spectroscopy extended the accretion tracers to the 
Pfund series lines, including the Pf/3 hydrogen recombination 
line at 4.65/rm (H 1(7-5)). 

To test if the H I mid-IR lines can be used as accretion in¬ 
dicators, in the rest of this section we focus on the subsample 
of full and transitional disks for which an estimate of the Ha 
equivalent width is available in the literature (see Table I Aik 
We chose to use the H 1(7-6) line (the strongest among the H I 
lines observed with the Spitzer SH module) to maximize the 
number of detected sources. Among the objects considered 





lo g L acc ( L o) 


Fig. 8 .— H 1(7-6) line luminosity versus Ha line luminosity and 
L acc . Circles are drawn around water-contaminated objects. The 
black-solid line shows the best linear regression line and the shaded- 
gray area represents the approximate 95% confidence intervals (2cr) 
on the regression line. The typical error is shown in the bottom-right 
comer. 


in this work, 73 have an estimate of the accretion luminosity 
from the Ha line (34 of those have an H 1(7-6) detection). 

We show in Fig. [8]the H 1(7-6) line luminosity as a function 
of the Ha line luminosity and t he accretion luminosity mea¬ 
sured as described in Sect. 13.21 There is a linear correlation 
between logL///( 7 _ 6 )-logL acc . To estimate the best fit regres¬ 
sion line for the objects with detection in both quantities, and 
to account for the measured errors in both the x- and y-axes 
we employe the Bayesian method as implemented in the IDL 
routine linmix_err (Kelly 2007). The error on the y-axes is 
driven by the uncertainty on the observed line fluxes. The un¬ 
certainty in our estimate of L acc is given by the combination 
of several factors: the uncertainty on the Ha equivalent width, 
the non-simultaneous observation of Ha, photometry, and ex¬ 
tinction, the error on the extinction value (A v ) and the errors 
on slope and intercept of the L acc -LH a relationship. Follow¬ 
ing Costigan et al. (2012) we assume a mean amplitude of 
variation in the derived accretion luminosities of ~0.37 dex, 
which includes Ha equivalent width uncertainties, time vari¬ 
ability, and the error in the L acc -LH a relationship. We also 
consider an uncertainty on Ay of 1 mag (e.g. Edwards et al. 
2013), which gives a final uncertainty on L acc ~0.5 dex. We 
find the following best fit regression line between the H 1(7-6) 
line luminosity and the accretion luminosity: 

logL HI a-6)/L e = (0.48 ±0.09) x logL acc /L e -( 4.68 ±0.10) 

(1) 

We have checked that the logLHi(i-6)-\ogL acc correlation is 
not driven by the distance and/or stellar luminosity depen¬ 
dence. We have also ascertained that the relationship between 
the luminosity of the H 1(7-6) line and L acc measured through 
UV-excess emission for only 15 objects (see Sect. 13.2b lies 
within the 95% confidence intervals. We note that the un¬ 
certainty on the intercept we find in Eq. (1) is similar to the 
those found using different H I lines (e.g., Alcala et al. 2014), 
suggesting that the scatter around the fit is similar when us¬ 
ing different H I lines. On the other hand, the slope found in 
Eq. (1) is shallower than the ones found for the H I lines in the 
optical and near-IR wavelengths. This suggests that the MIR 
lines are less sensitive to accretion luminosity variations. 

Knowing L acc and the stellar properties we can also com¬ 
pute the mass accretion rate ( M acc ) using the classical relation 















10 


Rigliaco et al. 



EW PV (A) 


Fig. 9.— Top: H 1(7-6) versus Pf /3 line luminosity. The black- 
dashed line hows the least square fit. Bottom: comparison between 
the equivalent widths of the two lines. 


from Gullbring et al. (1998): 





Rstar \ I L acc R star 

R in ' GM s t ar 


( 2 ) 


where R,„ is the radius at which the accreting gaseous disk is 
truncated due to the stellar magnetosphere. In full disks, R,„ 
has been found ranging between 3R star and 10R„ U ,. (Johns- 
Krull 2007), and usually R„,=5R Jfflr is assumed. In our sample 
we do not only consider full disks, but also transitional disks 
where the geometry of the gas disk around the star could be 
different. We do not speculate on how that could change for 
more evolved disks, and we will use the R,„=5R v/m - for all the 
different disk morphologies. 

For the subsample of objects where the Pf/3 line is detected 
(Salyk et al. 2013) and the H 1(7-6) is detected and not wa¬ 
ter contaminated (this work) we find that their luminosities, 
as well as equivalent widths, are correlated with Pearson cor¬ 
relation coefficients —0.7 and —0.8, respectively, and p-value 
—0.007 and —0.05, respectively (Fig. 0. While Salyk et al. 
2013 report lower Pf/3 luminosities and high EWs for 5 tran¬ 
sitional disks with respect to full disks, we do not find any 
statistically significant difference in the H 1(7-6) lines in our 
muc h lar ger samples of transitional and full disks (see also 
Sect. 14. 11 1. We also note that these two samples have not sta¬ 
tistically different L acc , as also reported in other studies where 
L flcc are consistently derived (Keane et al. 2014, Fang et al. 
2009). 

The fact that the H 1(7-6) line correlates both with the Ha 
and Pf/3 lines suggests that the Humphreys a line at 12.37/im 
is likely an accretion indicator, as the hydrogen lines belong¬ 
ing to other series (Hct (e.g., Muzerolle et al. 1998, Pa/3 (e.g., 
Folha & Emerson 2001), Br 7 (e.g., Natta et al. 2006), Pf/3 
(Salyk etal. 2013)). 

Speculating on the correlation reported in Eq. (1), we can 
measure what the mass accretion rates would be if the H 1(7- 
6) lines in the debris disks spectra were accretion related. We 
find that debris disks with detected H I lines span the lower 
range of mass accretion rate (covering the range between 
x 10 -10 -2 x 10 12 Mg/yr), as would be expected by their more 
evolved stage and the lack of gas in their disks. 


4.3. HI Line Ratios 


Hydrogen line ratios can be used as a diagnostic tool to con¬ 
strain the physical conditions (density and temperature) of the 
emitting gas. Among the sample of objects analyzed in this 
paper, in 46 the H 1(7-6) transition is detected above 3er, and 
11 of those show as well a detection of the H 1(9-7) transition. 
The detected H l(9-7)/H 1(7-6) line ratios range between —0.4 
and — 1.1. The upper limits cover a similar range. We can 
retrieve information on the physical properties of the emitting 
gas by comparing the observed ratios with those predicted by 
two different models: the Case B (Hummer & Storey 1987, 
Storey & Hummer 1995) and the Kwan and Fischer (Kwan & 
Fischer 2011, KF) models. In these models the H I line ratios 
ranges between —0.3-0.5 in the optically thin regime, and be¬ 
come higher (between — 1.0-2.0) when the lines are optically 
thick. 


4.3.1. Case B models 

The Case B model for radiative ionization and recombi¬ 
nation (Baker & Menzel 1938) assumes that gas is optically 
thick to Lyman series photons and optically thin to photons as¬ 
sociated to all other transitions. This model has been largely 
used to retrieve the properties of gas accreting onto T Tauri 
stars. However, results found considering individual objects 
and different hydrogen lines indicate very different physical 
properties of the accreting gas (Bary et al. 2008, Rigliaco et 
al. 2009, and Edwards et al. 2013 for a review). Using the 
interactive online server 3 (Hummer & Storey 1987, Storey & 
Hummer 1995), we computed the predicted line ratios for a 
range of temperatures from 1000 K to 15,000 K and electron 
densities from 10 9 -10 12 cm -3 (see Fig. 0. 

The predicted H l(9-7)/H 1(7-6) line ratios are almost con¬ 
stant (between —0.3-0.4) for T>5000 K over the electron den¬ 
sity range considered here. As the temperature decreases un¬ 
der 5000 K, the predicted ratio increases from —0.3 to —0.7 
until electron densities n e < 10 10 cm -3 , and then decreases to 
—0.2 when n e becomes larger. Thus, a large range of tem¬ 
peratures, down to T—1000 K, is required to explain most of 
the observed ratios. Still, ratios as high as 1.0 cannot be ex¬ 
plained by the Case B models. In addition, as pointed out in 
Edwards et al. (2013), the necessity to maintain low line opti¬ 
cal depths in the Case B models puts some restrictions to the 
applicability of this model to regions with high densities such 
as accretion funnels in YSOs. 

4.3.2. Kwan & Fischer models 

Kwan & Fischer (2011) predicted hydrogen line ratios us¬ 
ing a different approach than the Case B, in that they consider 
individual line optical depths. The inputs of these calculations 
are the gas temperature, the hydrogen density, the ionization 
rate (which was not specified in the case B models) and the ve¬ 
locity gradient with respect to the emission length scale (see 
Kwan & Fischer 2011 and Edwards et al. 2013 for more de¬ 
tails). These parameters are chosen to replicate conditions in 
the region of winds and accretion flows in accreting YSOs. 

We show in Figure QT| the comparison between the KF 
model predictions for the H l(9-7)/H 1(7-6) line ratio and the 
observations 4 . The input gas temperatures span from 5000 K 
to 15,000 K, and given that the (n,,)/(n//) ratio is between 0.1 
and 1.0, the hydrogen densities (n H ) cover similar electron 
densities as those shown in Fig. [TO] for the Case B models. 

3 http://cdsarc.u-strasbg.fr/viz-bin/Cat7VI/64 

4 Hydrogen line ratios for the Kwan & Fischer models are publicly avail¬ 
able: http://iopscience.iop.Org/0004-637X/778/2/148/suppdata/data.tar.gz 






Probing stellar accretion with mid-infrared hydrogen lines 


11 



Fig. 10.— HI(9-7)/HI(7-6) line ratio versus the electron density. 
Case B models for different temperatures as labeled. The range of 
the observed ratios is shown along the left side, where the typical 
uncertainties on the observed ratios are also plotted. 


The ratios suggested by the KF models are almost constant 
for all the temperatures until n// ~ 10 1() cm" 1 . As the density 
increases, the line optical depth increases as well, first for the 
(7-6) transition, which explains the steep increase in line ratio. 
For (n w )> 10 1 'em" 3 both lines have large optical depths which 
explains the plateau. While the general behavior is similar for 
all temperatures, the specific transition density is ~ 10 10 cm -3 
for T> 10,000 K and ~ 5 x 10 10 cm 3 for lower temperatures. 
Thus, unlike the Case B models, the KF models can explain 
all the observed line ratios with a limited range of tempera¬ 
tures (5000 - 15,000 K) and relatively narrow range of gas 
densities ( n H ~ 10 10 - 10 n cm -3 ). 



Fig. 11.— HI(9-7)/HI(7-6) line ratios versus the hydrogen column 
density for the observed ratios. Kwan & Fischer models for different 
temperatures as labeled. Symbols as in Fig. 1101 


4.3.3. Comparison between model predictions and gas 
physical properties 

According to Martin (1996), who modeled the tempera¬ 
ture and density profiles for gas channeled in magnetospheric 


accretion funnels, the typical temperatures in these regions 
are —6000 K, and hydrogen densities are as high as 10 9 - 
10 10 cm -3 . Muzerolle et al. (2001), considering the Ha and 
H/3 lines, found that the temperature of magnetospherically 
accreting gas ranges between 6000 K - 12000 K for a range 
of accretion rates between 10“ 6 - 10 10 M 0 /yr. Recently, Ed¬ 
wards et al. (2013) examined near-IR hydrogen line ratios 
for a sample of T Tauri stars. They find that the hydrogen 
density (hh) in the line formation region is constrained within 
2x 10 lo -2x 10 n cm -3 , while the temperature is not well delin¬ 
eated. 

If the hydrogen mid-IR lines are tracing accretion at least 
in full and transitional disks, the observed H l(9-7)/H 1(7-6) 
line ratios cannot be explained by the gas properties predicted 
by the case B models. In fact, as we have seen is Sect. 14.3.11 
and Fig. [lO] the observed line ratios require gas colder than 
5000 K irrespective of the electron densities. Such physical 
properties, and in particular such low temperatures, are un¬ 
likely in the magnetospheric accretion columns, in which we 
expect temperature of at least a factor of two higher and a 
much narrower range in density (e.g., Muzerolle et al. 2001). 

On the other hand, the KF models suggest gas properties 
that better match the scenario in which the mid-IR hydro¬ 
gen lines are produced in the accretion funnels because they 
predict an hydrogen column density («#) between 10 10 -10 11 
cm -3 , and a temperature range between 5,000 and 15,000 K 
(see Fig.Qj}. The analysis of the H l(7-6)/Pf/3 line ratios for 
the subsample of objects where these two lines have been ob¬ 
served adds additional constraints. In all but one object the 
observed line ratio is consistent with the KF model predic¬ 
tions for iih between — 10 10 — ~ 10 11 cm -3 and T>5,000 K 
(see fig.f]~2l). These results not only suggest that in most cases 
H 1(7-6) and Pf/3 might have a common origin, but that this 
origin is from hot gas with relatively high hydrogen column 
density (as expected in the accretion funnels). The objects 
where both H 1(7-6) and Pf/3 are detected, but the H l(7-6)/Pf/3 
line ratio does not fit the Kwan & Fischer models are VW Cha 
and TW Hya, where the H l(7-6)/Pf/3 line ratios is —0.9. 



Fig. 12.— HI(7-6)/Pf/3 line ratios versus the hydrogen column den¬ 
sity. Kwan & Fischer models for different temperatures as labeled. 
Same symbols as Fig. 1101 


We note that similarly stringent constraints on the gas phys¬ 
ical properties cannot be placed when using shorter wave- 
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lengths hydrogen lines. In particular, we have analyzed model 
predictions (both case B and KF models) for the H l(7-6)/Ha 
and H l(7-6)/Pa/3 lines. We find that the observed ratios span 
a wider range of values, making the comparison with models 
inconclusive. This could be due to the different optical depths 
of the compared lines, which could affect the Ha line as well 
as other short wavelength transitions, or to time variability, 
which may be more pronounced at shorter wavelengths. In 
the case of the H l(7-6)/Ha line ratio we must also take into 
account that Balmer series lines in general, and Ha line in 
particular, have quite complicated profiles, with a variety of 
structures that include blueshifted absorption and line asym¬ 
metries (e.g., Basri & Batalha 1990). 

4.4. Chromosphere and High-Density Wind 

The analysis carried out so far has shown that the H I mid- 
IR lines are very unlikely disk tracers. Moreover, the corre¬ 
lation between the H 1(7-6) line and the accretion luminos¬ 
ity, and the H I line ratios support an origin of these lines in 
the accretion funnels for full and transitional disks. However, 
Hollenbach & Gorti (2009) discussed two other possibilities 
for the origin of H I mid-IR lines: an internal wind shock 
with high speed (v s >100 km/s) and close to the wind origin 
(< 1 AU), or the stellar chromosphere (see also Pascucci et al. 
2007). 

An internal shock origin, even if plausible for full disks, 
would be very difficult to reconcile with the detection of H I 
in the spectra of transitional and debris disks. Transition disks 
rarely show signs of outflows activity (see for example the op¬ 
tical forbidden line profiles in Hartigan et al. 1995), and de¬ 
bris disks have already dispersed most of their gas disk mass 
(e.g. Pascucci et al. 2006). While it is possible that mid-IR 
hydrogen lines trace different environments depending on the 
stellar/disk evolutionary stage, as shown for the [Ne II] line 
(see Sect. [[]), a wind shock origin cannot explain the detec¬ 
tions in transitional and debris disks. 

The mean chromospheric electron density in T Tauri stars 
is n e ~ 10 n cm -3 (Giampapa et al. 1 981, o r 10 11 < n H < 
10 12 cm -3 if 0.1<(n e )/(ntf)<l, see Sect. 14.3.2] ). According to 
the KF model at these high densities the H 1(7-6) and (9-7) 
transitions would be both optically thick, hence line ratios 
should be >1. Only few sources in our sample sport such 
large line ratios. Another way that we have employed to inves¬ 
tigate the chromospheric origin of the H I lines is to consider 
the effect of the chromospheric emission on the measured ac¬ 
cretion luminosity. Manara et al. (2013), quantified this ef¬ 
fect by analyzing a sample of high dispersion spectra of non¬ 
accreting low-mass T Tauri stars. Measuring the Ha emission 
in non-accreting stars, and converting this value into an accre¬ 
tion luminosity, they find that the chromospheric activity can 
significantly contribute to the measured accretion luminosity 
when the latter is smaller that 10~ 3 Lq, and that this contribu¬ 
tion is strongly dependent on the object’s effective tempera¬ 
ture (hence spectral type). They derive a correlation between 
T e ff and L acc ,noise (the latter representing the chromospheric 
contribution to the accretion luminosity). Unfortunately, their 
sample is limited to objects with spectral type in the range K7- 
M9.5, and none of our debris disks falls in this range. In the 
following we will focus on the 42 obj ects with T e ff <4100 K 
and with an estimate of L acc (see Tahle lA I b . In Fig.[l3]we plot 
L acc,noise obtained using the correlation provided by Manara et 
al. (2013, their Eq. 2) versus L acc as reported in Table lAll 13 
of these objects also shown a detection in the H 1(7-6) line 
(green-filled symbols). There is a small number of targets. 


where the ratio is not more than 5, for which the chromo¬ 
spheric contribution could substantially influence the estimate 
of the accretion luminosity. However, in the majority of ob¬ 
jects L acc is at least 5 times bigger than its chromospheric con¬ 
tribution, meaning that the latter is well within the estimated 
uncertainty on L acc . Moreover, the objects where the H 1(7-6) 
line is detected (identified as green-filled symbols in Fig.fl3li 
span all over the range of L acc . 
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Fig. 13. — Contribution to the accretion luminosity from chromo¬ 
spheric emission (L acc fL aC c noise), measured using the relation found 
by Manara et al. (2013), versus L acc for all the objects with 
Tetf <4100 K. Symbols are as labeled. Filled symbols refer to ob¬ 
jects where H 1(7-6) is detected. 


We can only test the chromospheric contribution over the 
accretion luminosity for —40% of the total sample, where 
T e ff <4100 K. For this subsample the previous analysis sug¬ 
gests that the H I lines are not contaminated by the chromo¬ 
spheric contribution at a significant level. We could not extend 
this analysis to objects with higher temperatures (earlier spec¬ 
tral types) and to debris disks. It is possible that the H I mid- 
IR lines are tracing accretion in less evolved objects, while 
their contribution is only chromospheric in later evolutionary 
stages. 

The way to discriminate between accretion funnels, chro¬ 
mosphere and high-density wind is through ground-based 
high-resolution (10-20 km/s) spectroscopy. In fact, for H I 
lines originating in accreting gas we will observe broad lines 
(FWHM— 100 km/s, Muzerolle et al. 2001) close to the stel¬ 
lar velocity. A chromospheric component will be narrower 
(—40 km/s) and centered at the stellar velocity (e.g. Herczeg 
et al. 2007). An internal shock origin will instead produce 
H I profiles that are as broad as 100 km/s and can be also 
blueshifted by the same amount. 

4.5. A closer look at the Debris Disks 

As already mentioned in the previous sections, this paper 
reports the first detection of H 1(7-6) and H 1(9-7) gas lines in 
debris disks (we also detect the [Ne II] gas line in the Spitzer 
spectrum of HD 35850, as will be discussed in Rigliaco et al., 
in preparation). We show in Fig. [14] a compilation of the 8 
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Wavelengths (p, m) 


Fig. 14.— Spitzer spectra of the 8 Debris Disks where the H 1(7-6) 
line is detected. 


debris disks line profiles where the H 1(7-6) line has been de¬ 
tected. In the following we will investigate the properties of 
all the debris disks we have re-reduced, either with detected 
H 1(7-6) lines or not (as reported in Tables IAll and lA2l >. All 
these debris disks are observed within the FEPS 5 Spitzer Sur¬ 
vey (Meyer et al. 2006), and we collect the quantities plotted 
in Figure[l5]and reported in Table|4]from the FEPS Database 6 . 
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Fig. 15. — Top-left panel: X-ray luminosity versus Age. Black 
squares and red-open circles refer to objects with detected and not- 
detected H 1(7-6) line, respectively. Top-right panel: Lx/L star versus 
Ltf/( 7 -«)/Lj,ai-. Bottom-left panel: chromospheric activity diagnostic 
R hk versus the H I line luminosity. Bottom-right panel: H 1(7-6) line 
luminosity versus Age. 


TABLE 4 

Debris Disks Properties 


Name 

Age 

(Myr) 

log Lx 
(erg/s) 

r hk 

Ref 

1RXS J121236.4-552037 

17 

30.32 


|6l 

PI 


1RXS J122233.4-533347 

17 

30.34 


16 

2 


1RXS J130153.7-530446 

17 

30.01 


16 

2 


HD 119269 

17 

30.60 


p 

2 


1RXS J133758.0-413448 

15 

30.40 


|6 

2 


1RXS J161458.0-275013 

5 

30.05 


J7 

2 


HD 142361 

5 

30.71 

-4.02 

[712 

51] 

RXJ1600.6-2159 

5 

30.31 


f7 

2| 

RXJ1612.6-1859 

5 

30.65 


7 

% 

HD 35850 

12 

30.34 

-4.22 

1 1 012191 

V343Nor 

12 

30.33 


I10|2I 

HD 19668 

257 


-4.38 

12191 

HD 25457 

112 

29.74 

-4.39 

[fRDl 

HD 37484 

50 

29.16 

-4.49 

[2|2|5 

i] 

HD 12039 

63 

29.56 

-4.16 

bbl5 

i| 

HD 202917 

32 

30.07 

-4.22 

229 


HD 141943 

16 

30.62 

-4.36 

224 


HD 377 

87 

29.52 

-4.36 

22 9 


HD 17925 

55 

29.07 

-4.31 

2 2 1 


HD 32297 






HD 61005 

135 


-4.26 

|2]4 


HD 69830 

1860 


-4.95 

09 


HD 134319 

186 


-4.35 

09 


HD 216803 

400 

28.33 

-4.27 

[2l2l4l 

AO Men 

12 

30.16 

-3.75 

11012131 


References. — (1) Baliunas et al. 1996; (2) FEPS 

Database; (3) Grey et al. 2006; (4) Flem-y et al. 1996; (5) Luh- 
man et al. 2005; (6) Mamajek et al. 2002; (7) Preibisch et al. 
2002; (8) White et al. 2007; (9) Wright et al. 2004; (10) Zuck- 
erman et al. 2001; 


We first investigate a chromospheric nature for the mid-IR 
hydrogen lines. As shown in Figure[l5](top-left panel), all the 
debris disks where the H 1(7-6) line is detected are younger 
and with higher X-ray luminosity. The latter is a chromo¬ 
spheric indicator of magnetic activity (e.g. Feigelson et al. 
2007), and its decay with age could let us conclude that the HI 
line in debris disks are chromospherical, and as the magnetic 
activity decreases (with age) they are not detected anymore 
in older disks. Moreover, the X-ray emission of the debris 
disks with detected H 1(7-6) lines almost reaches the satura¬ 
tion level for chromospheric emission (Ly/L s , u ,- ~ 10“ 3 , e.g., 
Patten & Simon 1996, Figure[l5] top-right panel). This satu¬ 
ration level indicates that the chromospheric radiation emitted 
from the star cannot exceed a given fraction of the total stellar 
flux, and this limit is almost reached in all the observed de¬ 
bris disks, favoring again a chromospheric origin of the H I 
lines. Another test to check whether or not H I lines in debris 
disks are from chromospheric activity is shown in the bottom- 
left panel of Fig. [15] where we plot the stellar chromospheri¬ 
cal activity indicator R HK (collected from White et al. 2007). 
This indicator measures the core emission in the Ca II H and 
K lines, and the larger is the value, the more chromospheri- 
cally active is the star (Noyes et al. 1984) 7 . For the debris 
disks where this value has been measured there is no signifi¬ 
cant difference in terms of range covered by R HK as a function 
of H I line luminosity, showing that stars with the same chro¬ 
mospheric activity (similar R HK value) may or may not have 
detected H I line. On the other hand, the object with strongest 


5 Formation and Evolution of Planetary Systems 

6 http://feps.as.arizona.edu/index.html 


7 we note the reader that the R HK index may itself vary with time along the 
magnetic cycle of the star. 
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/ 

H I line has the largest R HK value, which would support chro¬ 
mospheric origin in this case. However, the sample is limited 
to only three objects, and statistical studies are not possible. 

We also speculate on the possibility that, as for full and 
transitional disks, the detection of the mid-IR H I lines in the 
Spitzer spectra of debris disks might trace accretion at very 
low levels. In fact. Fig. [15] bottom-right panel shows that the 
H I line is detected in the youngest systems from our sample, 
for which low level of accretion may be still ongoing. While 
multiplicity studies are not complete for debris disks, it is in¬ 
triguing to note that the two debris disks with a confirmed 
companion close en ough to contribute to the observed Spitzer 
spectrum (see Table lAll in Appendix A), have mid-IR hydro¬ 
gen line detections. This raises yet another possibility that the 
H 1(7-6) line observed in debris disks is due to an accreting 
or chromospherically active companion, instead of the debris 
disk object itself. However, this remains a speculation at this 
point given the low number of detections. 

In summary, even if a chromospheric origin appears the 
most likely based on the tests carried out so far, accretion can¬ 
not be ruled out. Only higher spectral resolution will be able 
to pin down the origin of mid-IR hydrogen lines. 

5. SUMMARY AND OUTLOOK 

We have re-reduced archival Spitzer /IRS spectra of 114 ob¬ 
jects in different evolutionary stages, from full to transitional 
to debris disks, and we have investigated the origin of the H I 
mid-IR lines. Our main findings are as follows. 

1. After correcting for the disk water contamination, we 
detect H 1(7-6) line in 26 full disks, 12 transitional and 8 
debris disks. The fainter H 1(9-7) line is detected in 6 full, 3 
transitional and 2 debris disks. 

2. We find H 1(7-6) luminosities similar to the [Ne II] 
12.81 ^m luminosities. These large luminosities exclude that 
the H 1(7-6) line traces the same hot disk surface probed via 
the [Ne II] line. 

3. Case B models cannot reproduce the observed H 1(9- 

7)/H 1(7-6) ratios except for T<3000K. On the other hand, 
KF models reproduce the observed ratios for reasonable 
assumption of temperatures and hydrogen column densities 
(n//). Using the KF models the n H is fairly well constrained 
between 10 10 < < 10 n cm _3 but the temperature is not. 

4. For the subset of full and transitional disks we report a 
positive correlation between the H 1(7-6) line luminosity and 
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the accretion luminosity. This trend, together with the results 
shown in point 3 suggest an origin of these lines from gas 
accreting onto the star. 

5. A pure chromospheric origin of the mid-IR H I lines, 
can only be tested for sources with T c yy <4100 K (40% of 
the total sample). For this subset of objects we find that the 
majority of them are not contaminated by chromospheric 
emission at a significant level. 

6. We report the first detection of H I lines in Spitzer 
spectra of debris disks, which occurs only the in the youngest 
systems studied here (< 20 Myr). A chromospheric origin of 
these lines seems likely, but it is not possible to definitively 
rule out either a long-lasting accretion signature nature, or 
the presence of a non-resolved accreting companion. 

Ground-based high-resolution spectroscopy will enable 
testing this origin. Extending accretion indicators at mid-IR 
wavelengths may enable measuring mass accretion rates in 
less evolved, still embedded Class I objects (where the star 
is still surrounded by a remnant infall envelope and accretion 
disk). Still, care should be taken in converting the mid-IR 
H I line luminosities into mass accretion rate, because the en¬ 
velope surrounding Class I sources could absorb part of the 
emission. More interestingly and on a speculative base, if the 
mid-IR H I lines were confirmed as accretion indicators in 
debris disks, these lines may be used to extend mass accre¬ 
tion rate measurements to very low values. If the detected H I 
lines are indeed tracing accretion, the luminosities we mea¬ 
sure for the 8 debris disks correspond to mass accretion rates 
< 1 (L 1 °Mg/yr. Optical and UV accretion diagnostics are not 
sensitive to such low mass accretion rates for sun-like stars 
(e.g., Manara et al. 2013). 

Finally, extending accretion indicators to the mid-IR will 
enable simultaneous observations of the [Ne II] line (the only 
confirmed photoevaporative wind indicator) and the H I line 
(an accretion indicator), making possible a direct comparison 
of the two most efficient disk dispersal mechanisms, photoe¬ 
vaporation and viscous accretion (Alexander et al. 2014). 
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hydrogen line ratios. E.R. is supported by the NASA’s As¬ 
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APPENDIX 

In Figure lATI we show the H 1(7-6) line for all the full and transitional disks where the line is detected. The H 1(7-6) line profile 
for debris disks are shown in Figure[l4] In Table lAll we report the source properties, and Table I A2l report the mid-IR line fluxes, 
errors and continuum levels around the three line of interest. 
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Fig. Al.— IRS spectra around the HI(7-6) line for the objects where the line is detected 
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Fig. Al. —continued. 


TABLE A1 
Source Properties 


Name 

Association 

Distance 

(pc) 

SpTy A v T eff 

(K) 

(M©) 

L star 

(L©) 

R-star 

(R©) 

HaEW 

(A) 

L acc Reference 

(L©) 

Full disks 


AA Tau* 
BP Tau 
CW Tau* 
CX Tau 
CY Tau* 
DD Tau ™-* 
DE Tau 
DG Tau 
DK Tau" 1 ’* 
DH Tau* 
DL Tau* 
DN Tau* 
DO Tau* 
DP Tau" 1 ’* 
DQ Tau 11 ’* 
DR Tau* 
DS Tau* 
FM Tau* 
FN Tau 
FT Tau* 

FZ Tau* 
HK Tau"' 


Taurus 

Taurus 

Taurus 

Taurus 

Taurus 

Taurus 

Taurus 

Taurus 

Taurus 

Taurus 

Taurus 

Taurus 

Taurus 

Taurus 

Taurus 

Taurus 

Taurus 

Taurus 

Taurus 

Taurus 

Taurus 

Taurus 


140 

140 

140 

140 

140 

140 

140 

140 

140 

140 

140 

140 

140 

140 

140 

140 

140 

140 

140 

140 

140 

140 


K7 
K7 
K3 
M2.5 
Ml 
M3 
Ml 
K6 
K7 
M0 
K7 
M0 
M0 
M0.5 
M0 
K5 
K5 
M0 
M5 
M3 
MO 
M0.5 


0.9 

0.5 

2.2 

0.8 

0.9 

0.4 

1.6 

1.4 

1.4 

1.0 

2.0 

0.4 

2.2 

0.9 

0.5 

0.5 

0.9 

0.9 

1.4 

3.8 

2.7 

3.0 


4060 

4060 

4730 

3490 

3705 

3415 

3705 

4205 

4060 

3850 

4060 

3850 

3850 

3850 

3850 

4350 

4350 

3850 

3125 

3415 

3850 

3850 


0.73 

0.73 

1.15 
0.36 
0.47 
0.29 
0.46 
0.92 
0.74 
0.56 
0.75 
0.55 
0.56 
0.56 
0.56 

1.16 
1.10 
0.57 
0.22 
0.31 
0.57 
0.56 


1.38 

1.62 

0.76 

0.62 

0.71 

0.13 

1.47 

1.15 
1.32 
0.76 

1.16 
1.62 
1.05 
0.85 
0.93 
0.96 
0.68 
0.52 
0.79 
0.16 
0.50 
1.12 


2.3 

2.4 
1.3 
2.0 
1.9 
1.0 
2.8 
1.9 
2.2 
1.8 
2.2 

2.7 
2.2 
1.9 
2.0 
1.9 

1.3 

1.5 

2.8 
1.0 

1.4 
2.2 


37.1 

40.1 
134.6 

19.8 
69.5 

181.8 
54.0 
84.0 

19.4 

53.4 
105.0 

11.9 
100.0 

85.4 

112.9 
106.0 

58.7 

70.9 

24.7 

254.3 
200.0 

29.3 


0.07 

0.07 

0.30 

0.003 

0.06 

0.05 

0.23 

0.17 

0.02 

0.03 

0.62 

0.016 

0.19 

0.05 

0.04 

0.49 

0.16 

0.02 

0.007 

0.49 

0.26 

0.028 


IC348-2 

IC 348 

315 

A2 

3.8 

8970 

2.80 

57.1 

3.1 



18 11 

IC348-6 

IC 348 

315 

G3 

3.9 

5830 

2.40 

16.6 

4.0 



T8TTT 

IC348-31 

IC 348 

315 

G6 

8.6 

5700 

1.60 

5.00 

2.3 



18 11 

IC348-37 

IC 348 

315 

K6 

2.8 

4200 

0.90 

1.30 

2.2 



18 11 

IC348-55 

IC 348 

315 

M0.5 

8.5 

3850 

0.60 

1.00 

2.2 



T8TTT 

IC348-68 

IC 348 

315 

M3.5 

2.1 

3740 

0.30 

0.50 

2.0 



18 11 

04216+2603 

Taurus 

140 

M0 


3500 

0.28 

0.02 

0.3 



Tim 

04385+2550 

Taurus 

140 

M0 

9.3 

3850 

0.57 

3.67 

4.0 

20.0 

0.50 

14 21 


HN Tau" 1 ’* 

Taurus 

140 

K5 

1.2 

4350 

0.78 

0.22 

0.8 

163.0 

0.15 

14 49 

IP Tau 

Taurus 

140 

MO 

0.5 

3850 

0.59 

0.17 

0.9 

11.0 

0.005 

14 49 

RW Aur™'* 

Taurus 

140 

K3 

0.4 

4730 

1.48 

1.60 

1.8 

76.0 

0.924 

14|49 

SU Aur 

Taurus 

140 

G2 

0.9 

5860 

1.85 

8.90 

2.7 

3.5 

0.14 

I 14 I 9 I 

UY Aur™ 

Taurus 

140 

K7 

1.1 

3850 

0.60 

0.91 

2.1 

72.8 

1.02 11419145 

AS 353A"' 

Perseus 

150 

M3 

2.1 

3415 

0.32 

0.68 

2.2 

124.5 

1.60 

|38|9| 

LkHa326 

Perseus 

250 

MO 

3.1 

3800 

0.49 

0.05 

0.5 

52.7 

0.10 


2071-HI 3/ NGC2068 

NGC2068 

400 

K1 

2.9 

5080 

1.68 

2.63 

2.0 

36.7 

0.42 

|11| 

2071-SI/ NGC2068 

NGC2068 

400 

K4 

4.1 

4600 

1.60 

4.10 

3.1 

3.7 

0.005 

[T 

T] 

RX Jill 1.7-7620™ 

Cha 

160 

K1 

1.5 

5080 

1.78 

3.09 

2.2 

7.6 

0.04 

43 

48 

SX Cha™'* 

Cha 

160 

M0.5 

1.1 

3850 

0.57 

0.61 

1.6 

33.2 

0.02 

15 

16 

SZ50 

Cha 

160 

M3 

2.1 

3350 

0.90 

0.78 

2.5 

46.0 

0.015 

15 

20 

TWCha 

Cha 

160 

MO 

1.2 

3850 

0.56 

3.10 

3.7 

28.3 

0.04 

15 

16 

VW Cha™'* 

Cha 

160 

K5 

2.4 

4350 

1.14 

4.02 

3.3 

71.7 

0.81 

15 

42 

VZ Cha* 

Cha 

160 

K6 

0.5 

4205 

0.93 

0.84 

1.6 

71.4 

0.099 

15 

42 
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WX Cha'"’* 

Cha 

160 

K7 

2.1 

4060 

0.77 

0.77 

1.8 

65.5 

0.07 

XX Cha'"’* 

Cha 

160 

M2 

0.6 

3560 

0.37 

0.18 

1.1 

133.5 

0.03 

Hen-3-600A 

TWHya 

56 

M3 

0.0 

3415 

0.26 

0.06 

0.6 

21.8 

0.007 

TWA 10 

TWHya 

100 

M2 

0.0 

5500 

0.90 

0.60 

1.0 



1RXS J132207.2-693812 

LCC 

86 

K1 

0.6 

5035 

1.23 

1.14 

1.4 

19.5 

0.07 

EX Lup* 

Lupus 

190 

M0 

0.0 

3850 

0.55 

1.31 

2.4 

31.3 

0.06 

GQ Lup* 

Lupus 

190 

K7 

0.4 

4060 

0.76 

4.53 

4.1 

31.7 

0.31 

HT Lup'" 

Lupus 

190 

K3 

0.0 

4730 

1.85 

4.00 

2.8 

7.3 

0.11 

RU Lup* 

Lupus 

190 

K7 

0.7 

4060 

0.74 

1.28 

2.2 

136.0 

0.69 

1RXSJ161411.0-230536'” 

Upper Scorpius 

140 

K0 

1.1 

5250 

1.62 

2.95 

2.0 



HD 143006 

Upper Scorpius 

169 

G7 

0.6 

5800 

1.12 

3.30 

1.7 

10.7 

0.19 

DoAr 24E'”’* 

pOphiucus 

160 

G7 

6.1 

5630 

2.15 

10.75 

3.3 

7.1 

0.24 

f DoAr 25* 

pOphiucus 

160 

K5 

2.9 

4350 

1.11 

2.19 

2.3 

12.0 

0.05 

SR21'” 

pOphiucus 

125 

G2.5 

9.0 

5950 

2.67 

28.0 

4.7 



WaOph6* 

pOphiucus 

160 

K6 

1.0 

4205 

0.90 

1.46 

2.3 

35.0 

0.02 

Cedi 10-IRS4 

pOphiucus 

160 

K2 


4800 

0.85 

0.40 

0.9 



Cedi 10-IRS2 

pOphiucus 

160 

K3 

3.3 

4730 

2.23 

9.64 

4.4 



Haro 1-16* 

pOphiucus 

160 

K3 

1.2 

4730 

1.80 

3.39 

2.7 

59.0 

0.24 

RNO 90* 

pOphiucus 

160 

G5 

1.0 

5770 

1.28 

2.74 

1.6 

76.0 

0.47 

EC 82 

Serpens 

415 

K7 

1.0 

4060 

0.75 

3.21 

3.3 

8.0 

0.03 

RXJ 1842.9-3532 

CrA 

136 

K2 

1.1 

4900 

1.33 

1.29 

1.5 

27.2 

0.06 

RX J1852.3-3700 

CrA 

136 

K3 

1.0 

4730 

1.27 

1.00 

1.5 

27.4 

0.06 

V1331Cyg 

Cyn 

550 

G8 

0.0 

5500 

2.81 

21.0 

4.8 

48.4 

2.01 

Transitional disks 

IC348-21 

IC 348 

315 

K0 

4.7 

5250 

1.60 

3.80 

2.4 

6.4 

0.08 

IC348-67 

IC 348 

315 

M0.75 

2.0 

3720 

0.50 

0.50 

1.8 

32.5 

0.03 

IC348-72 

IC 348 

315 

M2.5 

3.0 

3580 

0.40 

0.70 

2.1 

2.0 

0.0007 

IC348-133 

IC 348 

315 

M5 

3.6 

3200 

0.20 

0.20 

1.5 



DM Tau 

Taurus 

140 

Ml 

0.0 

3705 

0.47 

0.37 

1.3 

138.7 

0.046 

FP Tau 

Taurus 

140 

M5 

0.3 

3125 

0.22 

0.62 

2.5 

38.5 

0.001 

GM Aur 

Taurus 

140 

K3 

1.2 

4730 

1.27 

1.00 

1.5 

96.5 

0.58 

LkCa 15 m 

Taurus 

140 

K5 

1.0 

4300 

1.06 

0.74 

1.5 

13.0 

0.027 

St 34"' 

Taurus 

140 

M3 

0.2 

3415 

0.24 

0.01 

0.3 

51.6 

0.016 

UX Tau A'” 

Taurus 

140 

K5 

0.3 

4350 

1.17 

0.91 

1.9 

9.5 

0.03 

V836Tau 

Taurus 

140 

K7 

0.9 

4060 

0.81 

0.34 

1.2 

55.0 

0.16 

LkHa 330 

Perseus 

250 

G3 

1.8 

5830 

1.25 

2.78 

1.5 

20.3 

0.59 

2071-H9/ NGC2068 

NGC2068 

400 

K2 

1.6 

5250 

1.41 

2.03 

1.7 

1.2 

0.004 

CHX 22 (T54)" 1 

Cha 

160 

G8 

1.2 

5520 

1.69 

4.63 

2.2 



CS Cha 

Cha 

160 

K6 

0.8 

4205 

0.90 

1.50 

2.3 

54.3 

0.39 

SZ Cha'" 

Cha 

160 

K0 

1.9 

5250 

1.40 

1.90 

1.7 



Sz 18 (T25) 

Cha 

160 

M2.5 

1.6 

3560 

0.38 

0.32 

1.4 

8.0 

0.003 

Sz 27 (T35) 

Cha 

160 

K8 

3.5 

4060 

0.80 

0.47 

1.2 

100.0 

0.08 

TWHya 

TWHya 

56 

M0 

1.0 

3850 

0.57 

0.64 

1.7 

194.0 

0.09 

IM Lup 

Lupus 

190 

M0 

0.7 

3850 

0.56 

3.52 

3.9 

4.7 

0.02 

DoAr 21 

pOphiucus 

160 

K1 

6.6 

5080 

2.20 

27.3 

6.7 

1.5 

0.12 

Debris disks 

1RXS J121236.4-552037 

LCC 

108 

K0 

0.2 

5150 

0.90 

0.50 

0.9 



1RXS J122233.4-533347 

LCC 

124 

GO 

0.3 

5860 

1.24 

2.23 

1.3 



1RXS J130153.7-530446 

LCC 

108 

K2 

0.2 

4700 

0.82 

0.30 

0.8 



HD 119269 

LCC 

85 

G4 

0.1 

5700 

1.07 

1.31 

1.2 



1RXS J133758.0-413448 

UCL 

98 

K0 

0.3 

5250 

0.90 

0.68 

0.9 



1RXS J161458.4-27501 

Upper Scorpius 

114 

G5 

0.6 

5250 

0.90 

0.70 

1.0 



HD 142361'” 

Upper Scorpius 

101 

G3 

0.5 

5830 

1.72 

7.12 

2.5 



RX J1600.6-2159 

Upper Scorpius 

160 

G9 

0.8 

5410 

1.66 

3.77 

2.1 



RX J1612.6-1859'” 

Upper Scorpius 

127 

K0 

1.5 

5250 

1.50 

2.40 

1.8 



HD 35850 

/3Pictoris 

26.5 

F7 

0.0 

6000 

1.20 

1.78 

1.2 



V343 Nor”' 

/3 Pic tori s 

38 

K0 

0.0 

5250 

0.90 

0.67 

0.9 



HD 19668 

AB Dor 

40 

G8/K0 

0.0 

5500 

0.90 

0.60 

0.8 



HD25457 

AB Dor 

19 

F7 

0.2 

6000 

1.20 

2.50 

1.5 



HD37484 

Columba 

60 

F3 

0.1 

6740 

1.30 

3.90 

1.5 



HD 12039"’ 

Tuc-Hor 

42 

G4 

0.2 

6000 

1.00 

1.00 

0.9 



HD202917 

Tuc-Hor 

46 

G7 

0.1 

5750 

0.90 

0.70 

0.8 



HD 141943 

Field 

67 

GO 

0.1 

6030 

1.55 

6.19 

2.2 



HD377 

Field 

39.6 

G2 

0.2 

6250 

1.00 

1.50 

1.0 



HD 17925 

Field 

10 

K1 

0.0 

5080 

0.92 

0.55 

0.9 



HD32297 

Field 

112 

A0 

0.2 

8000 

1.90 

6.20 

1.3 



HD61005 

Field 

34 

G8 

0.1 

5759 

0.90 

0.60 

0.8 



HD69830 

Field 

13 

G8 

0.0 

5500 

0.90 

0.60 

0.85 



HD134319 

Field 

44 

G5 

0.0 

5770 

1.00 

0.72 

0.8 



HD216803 

Field 

8 

K4 

0.0 

4600 

0.66 

0.19 

0.7 



AOMen 

Field 

39 

K3 

0.0 

4435 

0.70 

0.26 

0.9 




26l42 


43129 

43129 


15153 


43151 


14121 


4l2 


412 


14121 


26150 

64177 



36142 


33131 


5 

7 

5 

7 

5 

7 

5 

7 

5 

7 

13 

5 

5 

1 

13 

5 


19 


51 


3 

5jfe 


5 E 

s 


|5 

71 

5 

7 

5 

7 

5 

7 

5 

7 

5 

7 

5 

7 
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TABLE A1 
Source Properties 


REFERENCES. — (1) Appenzeller et al. 1983; (2) Biller et al. 2007; (3) Brown et al. 2009; (4) Calvet et al. 2002; (5) Carpenter & Stauffer 2003 (6) Casali 
& Eiroa 1996 (7) Chen et al. 2014 (8) Cieza et al. 2010; (9) Cohen & Kuhi 1979; (TO) Duchene et al. 1999; (11) Espaillat et al. 2012; (12) Evans et al. 2012; 
(13) Fernandez et al. 1995; (14) Furlan et al. 2006; (15) Furlan et al. 2011; (16) Guenther & Emerson 1997; (17) Herbig & Bell 1988; (18) Herbst, W. 2008; 
(19) Hillenbrand et al. 2008; (20) Hughes & Hartigan 1992 (21) Kenyon et al. 1998; (22) Kraus et al. 2011; (23) Kraus & Ireland 2012; (24) Lafreniere et al. 
2008; (25) Lafreniere et al. 2014; (26) Lahuis et al. 2007; (27) Leinert et al. 1993; (28) Levreault 1988; (29) Mamajek et al. 2004; (30) Manoj et al. 2011 (31) 
Merin et al. 2010; (32) Metchev & Hillenbrand 2009; (33) Muzerolle et al. 2010; (34) Nguyen et al. 2012; (35) Petrov et al. 2014 (36) Pinte et al. 2008; (37) 
Pontoppidan & Dullemond 2005; (38) Prato et al. 2003; (39) Ratzka et al. 2005; (40) Rebull et al. 2010; (41) Reipurth & Zinnecker 1993 (42) Reipurth et al. 
1996; (43) Silverstone et al. 2006; (44) Strom et al. 1989 (45) Torres et al. 2008; (46) Uchida et al. 2004 (47) Walter 1986; (48) Walter 1992; (49) White & 
Ghez 2001; (50) White & Hillenbrand 2005; (51) White et al. 2007; (52) White et al. in prep; (53) Wilking et al. 2005; (54) Woitas et al. 2001 (55) Zuckerman 
& Song 2004 

NOTE. — Objects marked with * are water-contam inate d. "-Spectrosco pic Binary treated a s sin gle star 03. S tars marked with m are bin ary or multiple 
syst ems with the following separation: AS 353A 5. 6" i|22), CHX 22 0.25" {24}, DD Tau 0.56"(49)rDK Tau 2.4" {22}, DoAr 24E 2.05" {39V DP Tau 0.11" 
1221 . HD 12039 0.15" { 3 - (t his multiplicity has not been confirm ed <12I ). HD 142361 0.73" 132} - ( this multiplicity has not be en co nfirmed J32ll25t ). HK Tau 
2.4" {27). HN Tau 3.1" {27V H T Lup 2.8" {gV LkCal5 0.08 " {23V RX J1111.7-7620 0.27" {24k RX J1612.6-1859 0.21" {53, 1RXS J161411.0-230536 
0.22" {32}, RW Aur 1.4 " 1221. SR21 6.4 " 7TlV St 34 1.18" {52V SX Cha 2.2" {24}, SZ Cha 5.12" {24}, UX TauA 2.7,0.14" {Jo}. UY Aur 0.89" {27V 
V343NOR 10" {43, VW Cha 0.64,0.10 " {24pWX Cha 0.79 "{24V 
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TABLE A2 

Mid-IR line fluxes, errors and continuum level around the lines 


Name 

Flux a 

[Ne 11] 
lcr a 

cont b 

H 1(7-6) 

Flux a lcr a 

cont b 

H 1(9-7) 

Flux a lcr a 

cont b 


Full disks 

IC348-2 

7.66E-15 

1.81E-15 

0.24 

4.05E-15 

7.45E-16 

0.25 

<6.46E-15 


0.33 

IC348-6 

<2.23E-15 


0.078 

<2.73E-15 


0.077 

C1.18E-14 


0.11 

IC348-31 

<2.68E-15 


0.076 

<3.43E-15 


0.067 

<3.62E~15 


0.067 

IC348-37 

<2.46E-15 


0.055 

<2.22E-15 


0.059 

<2.92E-15 


0.081 

IC348-55 

2.16E-15 

7.28E-16 

0.034 

<1.94E-15 


0.037 

<3.31E-15 


0.052 

IC348-68 

<1.41E-15 


0.011 

2.68E-15 

4.97E-16 

0.011 

<3.51E-15 


0.012 

04216+2603 

<3.41E-15 


0.20 

<4.77E-15 


0.20 

<7.17E-15 


0.23 

04385+2550* 

4.33E-15 

1.66E-15 

0.35 

<7.61E-15 


0.34 

C1.08E-14 


0.38 

AA Tau* 

1.32E-14 

2.14E-15 

0.31 

<7.61E-15 


0.32 

<1.63E-14 


0.41 

BP Tau* 

<5.61E-15 


0.36 

<6.70E-15 


0.38 

<1.80E-14 


0.51 

CW Tau* 

<1.71E-14 


1.13 

3.45E-14 (5.73E-14) 

4.20E-15 

1.18 

<3.73E-14 


1.33 

CX Tau 

<2.38E-15 


0.17 

<3.01E-15 


0.17 

<5.55E-15 


0.20 

CY Tau* 

2.22E-15 

9.68E-16 

0.15 

<4.93E-15 


0.16 

<3.53E-15 


0.18 

DD Tau* 

2.02E-14 

4.73E-15 

0.93 

<1.37E-14 


0.96 

<5.02E-14 


1.14 

DE Tau 

<5.75E-15 


0.35 

7.81E-15 

1.40E-15 

0.35 

C1.34E-14 


0.49 

DG Tau 

2.43E-13 

1.13E-14 

4.64 

4.49E-14 

6.49E-15 

4.47 

<5.44E-14 


4.52 

DK Tau* 

<9.78E-15 


0.78 

<2.46E-14 


0.88 

<1.54E-14 


1.24 

DH Tau* 

<3.99E-15 


0.12 

<4.01E-15 


0.12 

<1.04E-14 


0.16 

DL Tau* 

1.24E-14 

2.15E-15 

0.87 

7.12E-15 (3.24E-14) 

2.09E-15 

0.88 

1.58E-14 (1.84E-14) 

2.52E-15 

0.93 

DN Tau* 

3.83E-15 

5.38E-16 

0.28 

<3.22E-15 


0.29 

<2.92E-15 


0.31 

DO Tau* 

<1.51E-14 


1.89 

<1.75E-14 


1.94 

<1.76E-14 


2.17 

DP Tau* 

1.16E-13 

3.63E-15 

0.63 

<2.42E-14 


0.66 

<2.44E-14 


0.78 

DQ Tau* 

9.06E-15 

3.33E-15 

0.60 

<1.14E-14 


0.61 

<1.19E-14 


0.64 

DR Tau* 

<2.12E-14 


1.88 

3.46E-14 (1.03E-13) 

6.93E-15 

1.88 

1.69E-14 (4.08E-14) 

8.65E-15 

2.10 

DS Tau* 

<6.33E-15 


0.24 

4.23E-15 (1.28E-14) 

1.10E-15 

0.25 

<6.32E-15 


0.36 

FM Tau* 

<2.26E-15 


0.25 

<4.85E-15 


0.27 

<6.71E-15 


0.34 

FN Tau 

5.44E-15 

1.23E-15 

0.68 

<5.71E-15 


0.67 

<5.59E-15 


0.83 

FT Tau* 

5.97E-15 

1.45E-15 

0.25 

<4.54E-15 


0.26 

<1.07E-14 


0.31 

FZ Tau* 

C1.75E-14 


0.98 

2.46E-14 (8.65E-14) 

3.42E-15 

1.02 

2.23E-14 (3.17E-14) 

7.12E-15 

1.15 

HKTau 

5.51E-15 

1.08E-15 

0.20 

3.44E-15 

7.21E-16 

0.20 

<4.42E-15 


0.24 

HN Tau* 

2.59E-14 

4.81E-15 

0.95 

<1.66E-14 


1.01 

<4.05E-14 


1.36 

IP Tau 

<6.22E-16 


0.013 

<8.16E-16 


0.014 

<1.94E-15 


0.020 

RW Aur* 

<2.00E-14 


1.43 

1.96E-14 (7.16E-14) 

5.92E-15 

1.46 

<3.23E-14 


1.76 

SU Aur 

1.78E-14 

3.63E-15 

2.14 

3.27E-14 

3.87E-15 

2.17 

<7.15E-14 


3.02 

UYAur 

<2.49E-14 


2.36 

<3.39E-14 


2.38 

<8.22E-14 


2.67 

AS 353A 

1.1IE-14 

2.50E-15 

1.11 

1.96E-14 

3.07E-15 

1.13 

C1.80E-14 


1.26 

LkHa326 

4.74E-15 

1.23E-15 

0.33 

<8.34E-15 


0.33 

C1.27E-14 


0.36 

2071-HI 3/ NGC2068 

1.01E-15 

2.43E-16 

0.024 

<9.83E-16 


0.026 

<1.71E-15 


0.038 

NGC20682071-S1 

<1.97E-15 


0.046 

<1.77E-15 


0.050 

<2.23E-15 


0.061 

RX Jill 1.7-7620 

9.97E-15 

1.28E-15 

0.22 

<5.96E-15 


0.23 

<9.55E-15 


0.27 

SX Cha* 

<7.08E-15 


0.49 

<1.34E-14 


0.51 

<1.42E-14 


0.65 

SZ50 

<3.20E-15 


0.20 

<4.78E-15 


0.21 

<7.16E-15 


0.22 

TWCha 

<4.35E-15 


0.062 

<4.01E-15 


0.055 

<6.63E-15 


0.042 

VW Cha* 

3.13E-14 

3.50E-15 

0.73 

3.74E-14 (5.98E-14) 

2.79E-15 

0.76 

1.91E-14 (2.64E-14) 

5.18E-15 

0.95 

VZ Cha* 

<5.92E-15 


0.37 

<8.84E-15 


0.40 

<8.61E-15 


0.46 

WX Cha* 

<5.08E-15 


0.28 

1.29E-14 (2.62E-14) 

1.34E-15 

0.30 

<1.01E-14 


0.43 

XX Cha* 

4.57E-15 

7.48E-16 

0.13 

<5.66E-15 


0.14 

<6.46E-15 


0.16 

Hen-3-600A 

<6.62E-15 


0.099 

<4.78E-15 


0.095 

<1.08E-14 


0.11 

TWA 10 

<2.44E-15 


0.018 

<2.09E-15 


0.019 

<3.05E-15 


0.022 

RX 132207.2-693812 

<2.59E-14 


0.75 

3.46E-14 

6.13E-15 

0.74 

3.17E-14 

1.23E-14 

0.84 

EX Lup* 

<1.63E-14 


3.13 

1.71E-14 (5.34E-14) 

6.85E-15 

3.21 

<4.43E-14 


4.12 

GQ Lup* 

<6.38E-15 


0.49 

4.39E-15 (1.23E-14) 

2.29E-15 

0.52 

<9.91E-15 


0.65 

HTLup 

<1.90E-14 


2.17 

<1.96E-14 


2.21 

<5.23E-14 


2.64 

RU Lup* 

4.54E-14 

8.41E-15 

4.29 

3.61E-14 (7.74E-14) 

8.24E-15 

4.36 

<6.13E-14 


4.80 

RX J161411.0-230536 

6.44E-15 

1.36E-15 

0.34 

<5.69E-15 


0.35 

<7.30E-15 


0.36 

HD143006 

<1.08E-14 


0.69 

<1.63E-14 


0.74 

<3.68E-14 


1.12 

DoAr 24E* 

1.38E-14 

4.20E-15 

2.54 

1.83E-14 (3.96E-14) 

7.58E-15 

2.58 

<4.14E-14 


2.96 

DoAr 25* 

2.87E-15 

9.27E-16 

0.21 

<7.78E-15 


0.22 

<3.84E-15 


0.24 

SR21 

<2.84E-14 


0.79 

<3.63E-14 


0.75 

<6.01E-14 


0.91 

Wa Oph6* 

8.55E-15 

1.83E-15 

0.82 

<2.25E-14 


0.82 

C1.65E-14 


0.92 

Cedi 10IRS2 

<4.43E-15 


0.11 

1.92E-15 

6.24E-16 

0.12 

<6.97E-15 


0.15 

Cedi 10-IRS4 

2.43E-14 

8.57E-16 

0.080 

<1.84E-15 


0.070 

<1.75E-15 


0.054 

Haro 1-16* 

<9.73E-15 


0.50 

1.05E-14 (2.16E-14) 

3.42E-15 

0.54 

<2.43E-14 


1.05 

RNO 90* 

2.10E-14 

1.08E-14 

2.05 

5.64E-14 (1.28E-13) 

6.22E-15 

2.10 

3.87E-14 (4.92E-14) 

9.83E-15 

2.49 

EC 82 

1.39E-14 

3.05E-15 

1.32 

<1.48E-14 


1.45 

<3.93E-14 


2.57 

RX J1842.9-3532 

4.75E-15 

7.21E-16 

0.092 

3.83E-15 

7.96E-16 

0.10 

<5.25E-15 


0.19 

RX J1852.3-3700 

8.66E-15 

4.28E-16 

0.035 

4.62E-15 

3.26E-16 

0.037 

C1.74E-15 


0.055 

V1331Cyg 

<3.30E-15 


1.01 

1.16E-14 

6.24E-16 

0.96 

<5.64E-15 


0.95 

Transitional disks 

IC348-21 1 

3.12E-15 

5.13E-16 

BB 

3.71E-15 

9.77E-16 

0.081 

<5.67e-15 


0.11 

IC348-67 

2.40E-15 

6.69E-16 

0.018 

<2.73E-15 


0.019 

<3.43E-15 


0.023 

IC348-72 

<2.93E-15 


0.011 

<2.04E-15 


0.011 

<2.12E-15 


0.012 

IC348-133 

<2.90E-15 


0.017 

<2.69E-15 


0.016 

<4.85E-15 


0.023 
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TABLE A2 

Mid-IR line fluxes, errors and continuum level around the lines 


DM Tau 

5.45E-15 

4.59E-16 

0.047 

<2.13E-15 


0.048 

<1.56E-15 


0.055 

FP Tau 

<4.73E-15 


0.067 

<3.76E-15 


0.070 

<4.31E-15 


0.078 

GM Aur 

1.02E-14 

7.11E-16 

0.097 

1.89E-14 

9.28E-16 

0.10 

1.20E-14 

1.02E-15 

0.17 

LkCa 15 

3.41E-15 

5.54E-16 

0.14 

2.94E-15 

6.71E-16 

0.15 

<7.89E-15 


0.27 

St 34 

1.40E-15 

5.66E-16 

0.029 

1.59E-15 

4.68E-16 

0.028 

<2.12E-15 


0.029 

UX TauA 

9.04E-15 

9.85E-16 

0.11 

9.02E-15 

8.53E-16 

0.11 

<4.98E-15 


0.14 

V836Tau 

<2.05E-15 


0.10 

<3.61E-15 


0.11 

<3.33E-15 


0.17 

LkHa 330 

<4.99E-15 


0.49 

8.84E-15 

1.91E-15 

0.50 

<2.49E-14 


0.69 

2071-H9/ NGC2068 

1.21E-15 

3.75E-16 

0.023 

2.44E-15 

3.84E-16 

0.024 

<2.06E-15 


0.038 

CHX 22 (T54) 

5.48E-15 

9.30E-16 

0.032 

<2.57E-15 


0.036 

<1.04E-15 


0.049 

CS Cha 

3.61E-14 

6.12E-16 

0.046 

6 .1 IE-15 

2.53E-16 

0.048 

3.07E-15 

7.43E-16 

0.064 

SZ Cha 

1.48E-14 

1.02E-15 

0.16 

6.10E-15 

1.32E-15 

0.17 

<6.93E-15 


0.25 

SZ 18 

7.63E-16 

2.33E-16 

0.0075 

8.68E-16 

2.14E-16 

0.0073 

<2.63E-15 


0.0080 

Sz 27 (T35) 

6.32E-15 

9.03E-16 

0.021 

3.03E-15 

7.27E-16 

0.022 

<2.67E-15 


0.026 

TW Hya 

5.01E-14 

2.07E-15 

0.61 

7.93E-14 

2.00E-15 

0.62 

3.34E-14 

3.64E-15 

0.81 

IM Lup 

7.03E-15 

1.54E-15 

0.48 

<4.58E-15 


0.49 

<6.37E-15 


0.60 

DoAr 21 

7.80E-14 

2.58E-15 

0.65 

<1.32E-14 


0.63 

<4.66E-14 


0.86 

Debris disks 

1RXS J121236.4-55203 

<6.09E-16 


0.014 

<8.71E-16 


0.016 

<8.1 IE-16 



1RXS J122233.4-533347 

<4.75E-16 


0.020 

8.56E-16 

1.95E-16 

0.021 

<7.81E-16 


0.026 

1RXS J130153.7-53044 

<4.95E-16 


0.0084 

<4.67E-16 


0.0092 

<8.19E-16 


0.011 

HD 119269 

<1.15E-15 


0.024 

<9.33E-16 


0.026 

<1.49E-15 


0.031 

1RXS J133758.0-413448 

<6.79E-16 


0.019 

5.44E-16 

2.04E-16 

0.021 

<8.21E-16 


0.025 

1RXS J161458.4-27501 

<8.43E-16 


0.0088 

<6.31E-16 


0.0098 

<8.00E-16 


0.012 

HD 142361 

<1.59E-15 


0.040 

2.98E-15 

1.47E-16 

0.043 

1.44E-15 

2.60E-16 

0.052 

RX 1600.0-2159 

<7.96E-16 


0.012 

1.20E-15 

1.91E-16 

0.013 

<7.63E-16 


0.015 

RX 1612.2-1859 

<6.28E-16 


0.028 

9.32E-16 

1.53E-16 

0.031 

<8.60E-16 


0.036 

HD 35850 

3.38E-15 

1.01E-15 

0.26 

9.02E-15 

9.27E-16 

0.27 

<6.30E-15 


0.34 

V343Nor 

<3.30E-15 


0.11 

4.81E-15 

8.09E-16 

0.12 

<5.64E-15 


0.15 

HD19668 

<1.14E-15 


0.051 

<1.62E-15 


0.056 

<2.56E-15 


0.067 

HD25457 

<1.86E-14 


0.53 

<1.53E-14 


0.58 

<1.86E-14 


0.68 

HD37484 

<2.66E-15 


0.080 

<2.16E-15 


0.086 

<4.47E-15 


0.10 

HD 12039 

<7.55E-16 


0.065 

<1.31E-15 


0.071 

<1.43E-15 


0.082 

HD202917 

<1.94E-15 


0.044 

<3.19E-15 


0.048 

<4.31E-15 


0.058 

HD 141943 

<1.27E-15 


0.074 

3.10E-15 

4.13E-16 

0.080 

2.97E-15 

4.91E-16 

0.094 

HD377 

<1.79E-15 


0.088 

<1.31E-15 


0.095 

<2.73E-15 


0.11 

HD 17925 

<3.42E-15 


0.61 

<7.36E-15 


0.66 

<8.26E-15 


0.78 

HD32297 

<3.26E-15 


0.054 

<4.28E-15 


0.053 

<5.09E-15 


0.052 

HD61005 

<1.49E-15 


0.066 

<1.32E-15 


0.072 

<2.40E-15 


0.084 

HD69830 

<4.54E-15 


0.56 

<4.75E-15 


0.62 

<1.05E-14 


0.81 

HD134319 

<2.21E-15 


0.047 

<2.55E-15 


0.051 

<2.73E-15 


0.060 

HD216803 

<5.66E-15 


0.74 

<7.99E-15 


0.81 

<6.96E-15 


0.95 

AOMen 

<1.04E-15 


0.051 

<1.03E-15 


0.056 

<2.08E-15 


0.065 


Note. — “Fluxes and errors are in erg s 1 cnT 2 . For the detections we report the la-error, upper limits are at 3a. b Continua are in Jy. Objects marked with 
an * refer to water-contaminated disks, and the fluxes reported are already corrected for the water contribution. In parenthesis we report the observed H I fluxes 
before the water subtraction. 














